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ABSTRACT 

Electron  field  emission  from  semiconductors  is  investigated 
both  theoretically  and  experimentally.  The  theoretical  predictions  of 
the  general  Stratton  theory  are  calculated  specifically  for  silicon,  in 
the  [lOO],  [lio],  and  [ill]  directions.  A  method  of  simplifying  the 
calculation  of  the  energy  distribution  for  arbitrary  semiconductor  bands 
is  obtained,  utilizing  the  effective  mass  approximation. 

>  Experimental  field  emission  energy  distributions  (FEEDs)  are 
reported  for  both  n-  and  p-type  samples  of  low  resistivity.  The  experi¬ 
mental  distributions  are  characterized  by  a  high  intensity  single  peak, 
of  energy  0.4  eV  or  more  below  the  Fermi  level,  with  a  subsidiary  peak 
of  lower  intensity,  rising  from  just  below  the  Fermi  level.  The  larger 
peak  drops  in  energy  with  increasing  field.  Presented  data  demonstrates 
that  this  peak  lowering  is  not  attributable^*/  sample  resistance.  Obser¬ 
vation  of  the  subsidiary  peak  is  linked  to  either  low  sample  tempera¬ 
ture  or  low  doping,  implying  that  the  carrier  concentration  affects  its 
presence. 

Experimental  FEEDs  are  compared  to  those  expected  theoreti¬ 
cally.  It  is  concluded  that  they  are  not  similar.  Comparison  with 
photoemission  work  indicates  that  the  large  peak  is  due  to  a  band  of 
surface  acceptor  states.  The  subsidiary  peak  is  tentatively  ascribed  to 
conduction  band  electrons. 
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>_  Finally,  a  phenomenological  model  of  photo- field  emission 
(PFE)  is  proposed.  Based  upon  both  FEED  and  ITE  experiments,  this  model 
assumes  that  emission  occurs  primarily  from  surface  states.  A  second 
component  of  the  current  is  due  to  tunnelling  of  photogenerated 
electrons.  In  addition  to  photoconductivity,  a  self- regulating  break¬ 
down  mechanism  is  necessary  for  qualitative  agreement  with  experimental 
data.  One  such  mechanism,  avalanche,  is  investigated  for  the  dielectric 
emitter  model.  Qualitative  agreement  is  obtained  with  the  characteris¬ 
tic  non-linear  Fowler-Sordheim  behavior  observed  experimentally. 
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CHAPTER  I 


INTRODUCTION 


Field  emission  (FE)  from  metals  has  been  studied  since  the 
1920's  [26].  In  1928,  theorists  Fowler  and  Nordheim  successfully  pre¬ 
dicted  the  FE  current- voltage  relationship  [20].  Experimental  verifi¬ 
cation  of  this  theory  and  utilization  in  work  function  studies  of  adsorp¬ 
tion  followed  [26].  Field  emission  energy  distribution  work  (FEED) 
confirmed  the  "total  energy"  theory  of  Young,  in  1959  [57, 60].  These 
theoretical  foundations  were  based  upon  the  free  electron  theory  of 
metals. 

Serious  experimental  challenges,  most  notably  the  anomalous 
FEED  from  tungsten's  ( 1 00 )  plane  [53],  were  met  with  a  more  compre¬ 
hensive  theory.  This  general  theory,  due  to  Stratton  [51,52],  takes 
account  of  bulk  band  structure  and  its  effects  on  the  tunnelling 
probability  of  non- free  electron  bulk  electronic  states.  This  theory 
was  considered  to  adequately  explain  the  experimental  results  [23]. 

Due  to  the  overall  success  of  the  Fowler-Nordheim  (F-N)  and 
Stratton  theories  for  metals,  work  was  extended  to  field  emission  from 
semiconductors.  Strong  deviations  from  F-N  behavior  were  demonstrated 
in  the  early  research  on  CdS  [4].  These  were  attributed  to  photocon¬ 
ductivity,  which  could  be  adapted  easily  into  the  FH  theory.  Later, 
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however,  Busch  and  Fischer  showed  that  work  function  measurements  from 
silicon  showed  an  inverse  relationship  to  that  expected — n-type  emitters 
demonstrated  a  higher  barrier  than  those  of  p-type  [ll].  Further,  the 
double  peak  with  sharp  energy  gap,  predicted  by  the  Stratton  theory,  was 
not  observed.  Instead,  two  peaks  of  varying  separation  were  generally 
seen  [34,36].  Also,  wide  energy  distributions  were  observed  even  from 
very  heavily  n-doped  samples,  contrasting  greatly  to  those  seen  for 
metallic  emitters  [28,30].  In  addition,  it  was  demonstrated  that 
photoconductivity  could  not  account  for  the  entirety  of  deviations  in 
emission  from  high  resistivity  samples  [36]. 

This  thesis  presents  a  study  of  field  emission  from  semi¬ 
conductors,  silicon  in  particular.  It  is  based  on  an  attempt  to  syn¬ 
thesize  the  data  of  previous  researchers  with  new  experiments  of  FEED 
and  photon  stimulated  field  emission  (or  Photo  Field  Emission,  PFE,  in 
the  literature) .  New  analysis  of  this  evidence  is  compared  to  both  the 
theoretical  predictions  of  the  Stratton  theory,  and  proposed  surface 
state  models.  In  particular,  the  implications  of  pbotoemission  studies 
on  Si  are  compared  to  the  experimental  results  on  FEED  from  semicon¬ 
ductors. 

Using  the  results  of  this  comparison,  a  new  model  of  PFE  is 
developed.  Although  the  model  utilizes  crude  approximations  of  the 
quantities  involved,  the  predictions  of  the  model  show  some  qualitative 
agreement  to  the  actual  experimental  data.  This  model  is  meant  to  be  an 
incipient  attempt  to  link  together  the  experimental  parameters  of  the 
PFE  problem. 

In  undertaking  the  thesis,  novel  methods  of  sample  prepara- 
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tion  and  data  acquistion  were  developed.  Photolithographic  techniques 
were  applied  to  precisely  define  sample  dimensions.  A  suitable  pre¬ 
cleaning  treatment  was  found  to  reduce  surface  desorption  problems 
usually  associated  with  strong  oxide  layers  on  the  Si.  Also,  an  adapted 
van  Oostrum  type  electron  energy  analyser  was  interfaced  to  an  Apple  II 
microcomputer.  The  computer  control  allows  the  user  a  number  of  fea¬ 
tures,  such  as  rapid  data  acquistion,  graphing,  and  numerical  differen¬ 
tiation,  which  previously  were  performed  manually.  These  techniques  are 
described  in  detail. 

The  application  of  the  Stratton  theory  to  the  silicon  band 
structure  is  developed  in  an  original  way,  using  a  simplified  general 
approximation  to  the  constant  energy  shadows.  This  calculation  reduces 
the  Stratton  integrals  to  a  form  amenable  to  rapid  numerical  calcu¬ 
lation.  Direct  application  is  made  for  emission  from  the  silicon  (100), 
(110),  and  (ill)  planes. 

The  thesis  begins  with  a  synopsis  of  Young’s  free  electron 
FEED  theory.  This  second  chapter  also  contains  a  section  comprising  an 
original  extension  of  Stratton' s  general  formalism  to  cover  the 
particular  case  of  silicon.  This  provides  a  theoretical  basis  for 
comparison  of  the  experimental  FEEDs  which  will  follow. 

The  third  chapter  introduces  the  experimental  apparatus  and 
electronics  developed  for  FEED  studies.  Some  brief  discussion  of  the 
software,  also  originated  for  this  work,  follows.  Furtner,  improvements 
and  innovations  in  the  sample  preparation  are  discussed.  The  equipment 
used  in  PFE  studies  is  also  described  in  this  chapter. 

Experimental  results  are  explained  in  the  fourth  chapter.  The 
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first  section  of  this  chapter  summarizes  results  obtained  from  FEED, 

PFE,  and  F-N  techniques,  and  gives  surface  micrographs  for  samples  of 
various  dopings  and  temperatures.  Subsequent  sections  of  this  chapter 
compare  these  results  to  those  of  previous  researchers.  The  last 
section  discusses  a  generalization  of  these  observations  by  comparison 
to  the  photoemission  work  of  Wagner  and  Spicer. 

The  fifth  chapter  discusses  PFE.  Briefly  deriving  the 
photoconductivity  additions  to  the  F-N  equation,  it  is  demonstrated  that 
these  simple  theoretical  additions  cannot  completely  satisfy  the  experi¬ 
mental  data.  Using  the  conclusions  of  the  fourth  chapter,  a  new  model 
is  developed,  which  yields  more  satisfying  qualitative  behavior. 

The  sixth  chapter  comprises  a  general  conclusion  to  the 
thesis.  Here,  a  brief  summary  is  presented  of  the  previous  discussions. 
Following  the  conclusion  are  appendices  containing  the  computer  programs 
utilized  in  this  study. 
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CHAPTER  II 


THEORY 


Field  Emission  from  Metals 


The  general  theory  of  field  emission  from  free  electron  metals 
vas  developed  in  1928  by  Fowler  and  Nordheim  [20].  Their  solution  of  the 
current- field  relationship,  the  Fowler-Nordheim  (F-N)  equation,  provided 
one  of  the  first  applications  of  wave  mechanics  to  a  tunnelling  prob¬ 
lem  [26]. 

Young  rederived  the  F-N  equation  along  with  the  total  energy 
distribution  (TED)  in  a  classic  paper  [57].  The  formalism  is  repeated 
here,  as  silicon  energy  distributions  will  later  be  compared  to  the  ap¬ 
proximately  free  electron  distributions  from  tungsten. 

From  Figure  1 ,  we  see  that 


V(z)  -  -V 

z<0 

(1) 

a 

■  -ke2/4z  -  eFz 

z>0 

(2) 

where  the  first  term  of  (2)  is  due  to  the  image  potential  energy,  and 

the  second,  the  applied  field.  Here  we  distinguish  between  energy 

normal  to  the  surface  V,  and  total  energy,  E  by 

W  -  E  -  p^2/2m  -  Py2/2m 

-  p  2 /2m  +  V(z) 
z 


(3) 


TOTAL  ENERGY 


Figure  1.  Energy  Diagram  for  Metallic  Qnitter  with  Applied  Field 
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where  p.  is  the  momentum  in  the  ith  direction,  and  m  the  free  elec- 
tron  mass. 

By  defining  the  following  terms.  Young  placed  the  problem  on 
managable  ground: 

N(W,E)dWdE  “  number  of  e~  incident  at  the  surface  in  the 
interval  (E,E+dE) ,(W,W+dW)  per  unit  time  per  unit  area. 

D(W)  *  probability  of  barrier  penetration. 

P(V,E)dWdE  -  N(W,E)D(W)dVdE  =*  number  of  e~  in  interval 
which  penetrate  barrier. 

P(E)dE  *  (  /P(W,E)dW)dE  ■  Total  Energy  Distribution. 

j  *  e  /p(E)dE  *  current/ area  *  P-N  equation. 

Now  we  determine  these  quantities. 

From  Figure  2,  we  see  that  the  component  of  v  in  the  z  di¬ 
rection  is|"v|cos©.  The  relative  area  encompassed  goes  as  sin©,  so 

N(w,E)dwdE  -  n(E)dE[  |v|cos0sin0d9do/4rr]  (4) 

where  at  is  the  solid  angle,  4*  (steradians)  the  normalization,  n(E)  the 
density  of  electrons.  Then 

V  *  m(  Ivl  cos©)2)/2  (5) 

and 

-dW/m|v|  *  | v| sin©cos©d©  -  -dW/[2m(E-V)]1  ^  (6) 

After  substitution  into  (4)  and  integration  over  o, 

N(oi,E)dadE  -  -dWdEn(E)/2[2m(E-V)]1//2  (7) 

using  n(E)dE  *  4*(2ap/2 (E-V ^2dE 

h^(exp[ (E-u)/kT]  +0  (8) 

gives  N(W,E)dWdE  “  -4^m  dWdE 
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Using  a  tunnelling  probability  D(w)  of  the  form 

D(W)  -  exp[  -c  +  (¥  -w)/d]  (10) 

where  o  *  A(2m<p  J1  x 

v((ke5F  )^2/v)  (11) 

d  *  •fieP/(2(2m(9)1 /2  x 

t(  (ke^y)1  ^2/o)  (12) 

Both  t(y)  and  v(y)  are  tabulated  functions  of  value  near  1  ,  and  are 

discussed  following  equation  (20)  in  the  next  section. 

Upon  integration:  ^ 

_WA 

P(E)dE  -  ^N(W, E)B (W)  dWdE  (13) 

W=E 


■  4wmd  exp(-c-,u/d)  exp(E/d)dE 

- —  x  - 

h')  1  +  exp([E-i»]/kT) 

Then  the  total  current  density  j  is: 

» 


j  -  e  j P(E)dE 


(H) 


05) 


*  [e3F2/8:r  hot  ((a)1/2/p)]  x 

exp[[  -4(2m)1/%3,/2/(3fieF)]  x  [ v(  (a) 1/2 /»]]  x 
(*kT/d)/(sin(*lcT/d))  (16) 

where  a  *  (ke^y). 

Energy  distributions,  typical  of  the  free  electron  model  de¬ 
rived  above,  have  been  observed  by  many  researchers.  Various  distri¬ 
butions  appear  in  Figure  3,  as  predicted  by  Young's  equations  above. 
These  show  that  the  distribution  widens  with  increasing  temperature. 
Although  this  fact  has  been  borne  out  by  experimentation,  there  have 


I 
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Figure  3.  Theoretically  Predicted  Free  Electron  Total  Energy 
Distributions. 


dJ/dE 
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also  been  found  crystallographic  directions  of  tungsten  and  molybdenum 
which  do  not  exhibit  the  free  electron  distribution  predicted  [53].  To 
explain  the  phenomenon  of  "anomalous"  energy  distributions  from  metals, 
Stratton's  theory  of  field  emission  from  non- free  electron  solids  was 
invoked  [53].  This  theory  takes  account  of  the  actual  band  structure  of 
the  sample  in  describing  the  transverse  energy,  and  will  be  described  in 
the  next  section. 


Field  Emission  from  Semiconductors 


General  Formalism 

The  free  electron  calculation  of  Young  was  generalized  to  non- 
free  electron  solids  by  Stratton  [51 >32].  This  generalization  has 
apparently  been  quite  successful  in  explaining  the  anomalies  obtained  in 
TED's  from  tungsten  in  particular  emission  directions  [23>53]>  The 
Stratton  formalism  has  since  been  the  general  theory  for  field  emission 
from  bulk  states.  Although  measured  TEDs  from  semiconductors  do  not 
agree  with  this  theory  [5>30],  it  provides  a  theoretical  basis  for  the 
TED  expected  from  bulk  states.  For  this  reason,  the  present  section 
comprises  an  investigation  of  the  predictions  of  this  theory. 


Stratton  has  pointed  out  that  equation  ( 1  5 )  may  be  written 

[52]: 

P(E)  -  [2ef(E)/h3]  jf  D(V)  dpy  dpx  (17) 

or  equivalently: 

P(E)  -  [2ef(E)/«(2*)3]  jj  D(V)  dky  dkx  (18) 

where  the  integration  ranges  over  the  shadow,  on  the  emission  plane,  of 
the  constant  energy  surface  [23].  This  form  is  useful  for  semiconductor 
calculations  involving  subbands  or  multivalley  bands.  The  P(S) 
calculations  for  each  component  sub band  may  be  added  algebraically  [47]. 
To  proceed  further,  one  must  know  the  explicit  band  structure  of  the 
sample  in  question.  That  of  Si,  for  example,  is  shown  in  Figure  4. 

To  use  the  Stratton  formalism,  one  must  first  derive  expres¬ 
sions  for  the  transverse  energy  E  .  This  transverse  energy  is  unavail- 

£ 

able  for  tunnelling,  and  must  be  subtracted  from  the  total  energy  E  in 
the  integral  over  D(W) .  The  procedure  is  straightforward  in  the  <100> 
directions,  for  the  constant  energy  ellipsoids  of  the  conduction  band 
are  either  along  or  normal  to  these  axes.  For  other  directions,  the 
procedure  for  finding  E  is  non- trivial. 

w 

Further,  the  form  of  D(w)  given  in  equation  (10)  is  not  gener¬ 
al.  The  form  used  by  Young  was  an  approximation  of  the  general  VKB  re¬ 
sult,  with  image  corrections,  valid  near  E?.  The  general  result 
for  metals  is  [23]: 

D(V)  -  exp[-(4/3eF)(2m/fl2)l/2 

x  O  -  (V  -  EF)]3/2 

1/2 

x  vttke3?)  /  (®  -  (W  -  Ep)]  ]  (19) 

D(V)  -  exp[  -0.685  (®-  (W  -  Ef))5//2/F 


or 
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x  v(3.79  ?V2/  (®-  (W  -  BF)))  ]  (20) 

with  F  in  V/A,  o ,  W,  E  ,  in  eV.  In  the  above  equations,  v(y)  and 

r 

t(y)  can  be  approximated  by: 

v(y)  =  AO  +  Aly  +  A2y2  +  A3y5  (0<y<1  ) 

with  AO  »  1.0056,  A1  -  -0.1521,  A2  -  -1  .0625,  A3  =*  0.2083 
and  t(y)  5  1  +0. 1y 

The  correct  form  for  semiconductors  of  dielectric  constant  iC  differs 
slightly,  due  to  reduction  of  the  image  potential  by  the  factor  [(K-l)/(K 
+1  )]17^  ,  as  discussed  by  R.  Gomer  [25]. 


General  Band  Considerations 

Although  the  actual  integrations  must  be  performed  to  deter¬ 
mine  P(E)  for  any  particular  emission  direction,  the  transverse  energy 
at  the  conduction  band  edge  may  be  expressed  exactly.  For  arbitrary 
emission  direction  I)  *  hi  +  kj  +  lk,  there  exist  three  sets  of  2-fold 
degenerate  lobes — one  pair  per  axis.  The  centers  of  these  lobes  are  at 
an  energy  of  3.7eV  along  each  axis,  from  the  origin  [ll].  Thus  the 

transverse  energy  along  the  emission  direction  is  3*7eVsin  with  &  the 

.  2  2  2  2  2  2 
angle  between  D  and  the  axis.  But  sin  B  *  (k  +  1  )/(h  +  k  *  1  ) 

for  the  i  axis,  and  the  same  for  the  other  axes  by  cyclic  permutation  of 

the  indices.  Thus 

S  (hkl)  -  3.7eV  (k2  ♦  l2  )/(h2  +  k2  ♦  l2  )  (21) 

6  X 

E  .(hkl)  -  3.7eV  (h2  *  l2  )/(h2  ♦  k2  ♦  l2  ) 

*  j 


(22) 
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Figure  5.  General  Mapping  of  Constant  Energy  Ellipsoids  onto 

Emission  Plane.  Here,  a  ana  a'  denote  the  semimajor  axes, 
b  and  b'  denote  the  semiminor  axes.  Unprimed  variables 
correspond  to  the  constant  energy  ellipsoid,  primed  ones  to 
the  shade*  on  the  emission  plane. 
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E  (hkl)  -  3.7eV  (h2  +  k2  )/(h2  +  k2  +  l2  )  (23) 

b  fC 

The  general  solution  of  P(E)  for  the  silicon  conduction  band 
may  be  made  with  the  following  assumptions.  First,  because  of  the 
symmetry  of  the  constant  energy  ellipsoids  about  their  long  axis,  the 
angle  between  the  long  axis  and  the  emission  vector  D  completely  speci¬ 
fies  the  shape  of  the  transverse  energy  shadow.  The  dimensions  of  this 
shadow,  of  course,  must  come  from  E  and  the  effective  masses. 

Secondly,  it  is  assumed  that  each  3-d  constant  energy  ellip¬ 
soid  is  mapped  into  a  2-d  ellipse  on  the  surface  plane.  This  assumption 
allows  us  to  use  the  construction  in  Figure  5,  whereby  we  must  compute 
the  maximum  transverse  distance  in  the  2-d  k-space  of  the  emission  plane 
to  determine  a'  for  the  2-d  ellipse.  This  a'  in  turn  will  give  the 
"effective"  effective  mass  for  the  2-d  ellipse  on  the  long  axis.  The 
value  for  b'  of  the  2-d  ellipse  will  always  be  given  by  the  same  value, 
b.  This  is  intuitively  true  because  the  width  of  the  shadow  of  the  3-d 
ellipsoid  is  constant,  although  the  length  of  the  shadow  will  depend 
upon  the  angle  0. 

We  begin  by  describing  the  projection  of  the  ellipsoid  in  the 
D-i  plane  by 

K(kx,ky)  -  k„2/a2  +  ky*/b2  -  1  (24) 

with  a2  -  2m  E/-ft2  ;  b2  >  j/-fi2 

X  c 

putting  the  equation  of  the  ellipse  into  standard  form.  Here,  m^  and 
*  are  the  longitudinal  and  transverse  effective  masses,  of  value  0.91m 
and  0.19m  respectively  [8].  We  wish  to  maximize  the  distance  d(k„,ky)  be¬ 
tween  a  point  kx,ky  and  the  emission  vector  D  ■  hi  ♦  kj  ♦  lk,  given  by 
d(k„,ky)  -  k„sin&  -  kycos0  (26) 
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Using  the  Euler-Lagrange  multiplier  method,  we  take 


<Sd(kx,ky)  =*  A<5K(kx,ky) 

(27) 

Akxsin©  -  ikyCOs©  =  A(2kxAkx/a2  +  2kyAky/b2) 

or 

(28) 

(2kx.l/a2  -  sin©)Akx  +  (2kyA/b2  +  cos©)_lky  =  0 

(29) 

for  any  _lkx,  Aky.  Therefore,  each  expression  in  parentheses  in  (29) 

must  be  zero,  so 

kx*  *  a2sind/2A  ;  ky*  ■  -b2cos©/2A 

Substituting  (30)  into  (24)  yields 

kx*2/a2  +  ky*2/b2  *  1 

(30) 

so  A  »  ±(a2sin2©  +  b2cos2©)1/2/2 

(31) 

We  choose  the  +  sign,  since  maximal  distance  is 

kx>0,  ky<0.  Substituting  into  (30)  again,  gives 
kx*  -  a2sin©/( a2sin2©  b^os2©)1  ^2 

obtained  when 

ky*  --b2cos©/(a2sin2©  +  b^cos2©)1  ^2 

Again  substituting  (32)  into  (26)  yields 

(32) 

d(kx*,ky*)  *  (a2sin2©  +  b^os2©)1  ^2  ■  a’ 

The  new  2-d  ellipse  must  be: 

(33) 

k  2/a'2  ♦  k  2/b'2  -  1 
t  z 

with 

(34) 

a’2  -  2m*S/fi2  -  a2sin2©  +  b2cos2© 

(35) 

b’2  -  2m  S/1S2 

Installing  a  and  b  from  equation  (24)  yields  for  m*: 

m*(©)  *  m  sin2©  +  m  cos2©  (36) 

X  w 


The  general  expression  for  conduction  band  transverse  energy  must  take 
into  account  the  shifted  origin  of  the  2-d  ellipse.  As  described  in 
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equations  (21 )  to  (23)  this  implies  that 

dl2k  2  /(2m*(©))  -  3.7eV  sin2© 
to 

with  k  =  [(3.7eV)2m*(©)  sin2©  /  IT2]1 /2  (37) 

to 

and  a  transverse  energy  described  by 

E.(e)  =  H2(k.  +  k.  )2/  2  m*  (© ) 
t  t  to 

♦-aak2/(2mj  (38) 

Z  t 

Finally,  we  may  generalize  the  P(E)  from  the  conduction  band  to: 

P(E)  -  2  x  [  2ef(E)/«(2*)3  ]  x 

3  .  •  r  bl(kz) 

(  /  I  /  D(E  -B„  )  dkt|]  dk  \  (39) 

~b|(kz)  J  ’ 

a  -  [2m  E/  -H2]1 /2  (40) 

t 

ain2©,  -  [h2(1  -  5,, )  +  k2(l  -  <S,2)  ♦  1*0  -  ^f3)]  x 

1/(h2  +  k2  +  l2  )  (41a) 

cos2©,  -  (<S(1h2  ♦  <Si2k2  ♦  i|3la)/(h2*  k2  ♦  l2  )  (41b) 

<*(1  -  1  if  i-j;  0  if  i*j 

n,#  *  n^ain2©  ♦  m^cos2©  (42) 

kol  -  [2  (3  7eV)  m,*sin2©  /  dl2]1 /2  (43) 

b,  -  [^m,*/*2)  x  (E  -lAt*/20]1/2  (44) 

V 

3„  •  Mk„  ♦  kol  )2/(2fl»,  )  ♦  d^k  */ (2m  )  ]  (45) 

It  should  be  noted  here  that  the  shadow  of  each  constant  energy  ellipsoid 
is  taken  into  account  separately,  even  if  the  shadows  overlap  on  the 
emission  plane.  Calculations  for  the  three  major  crystallographic  planes 
are  described  in  particular  in  the  following  sections. 


! 
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Specific  Conduction  Band  Calculations 


lOO]  Direction 

The  general  form  given  above,  equations  (39)  to  (45),  can  be 
broken  into  two  parts  for  the  [lOO]  direction.  These  two  cases  arise 
because  of  the  degeneracy  in  energy  of  the  ellipsoids  symmetrically 
centered  on  the  j  and  k  axes,  as  shown  in  Figure  6,  and  the  ellipsoids 
doubly  degenerate  in  momentum  with  centers  on  the  i  axis.  It  can  be 
shown  that  (39)  simplifies  to 

P(B)  -  [2ef(E)/(dS  (2*)3)  j  x 


a  b 


(4 j  (y  D ( E  -  E,)  dkt)  dk*+ 

-a  b 

c 

2  x  (2t)  (  D(E  --ft2kt2/(2mt))  k,dk,  ) 

J 

(46) 

0 

c  -  [2m  E/ 

Um^E /  -K2  -  m1kl2/mt]1  /2 

(47) 

+  k0)2/(2ml)  ♦  -fi2k2/2mt 

(48a) 

[2  (3.7eV)m1/-h2]l/2 

(48b) 

[I10l  Direction 

The  [110]  direction  conduction  band  of  silicon  exhibits  2- 


fold  momentum  degeneracy  of  ellipses  centered  on  the  i  and  j  axes.  As 
shown  in  Figure  7,  ellipsoids  centered  on  the  positive  side  of  these 
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axes  lie  above  those  centered  on  the  negative  side.  These  overlapped 

ellipsoids  3till  contribute  to  the  transverse  energy  integrals  over  lc  , 

t 

It  ,  though  their  shadows  are  degenerate  with  those  of  the  ellipsoids 
z 

above.  For  this  reason,  the  form  of  P(E)  for  the  conduction  band  is: 
P(E)  -  [2ef(E;/«(2*)5]  x 

a  b1 

(4  J[  f  D(E  -  Et1  )  dkt1]  dkz  + 

-a  -b, 

i  bj 

2  J[  j  D(E  -  Et2)  dk,2]  dkx  )  (49) 

-a  -b2 

a  -  [2m  E/-K2]1  /2 

b 

b,  -  ([2m,*/*2]  X  [E  -  -fi2k ,2/2m^]  J1  (50) 

m,*  »  (mi  +  ®t)/2  ;  m2*  -  (51 ) 

k01  “  [■,*  3.7eV/*2]1/2 

koa  ’  [z^.w/n2]1^  (52) 

Et|  -  [-R 2( lc t ,  *  k0,  )2/2m,*  +  -h2k2/2mt]  (53) 

[ill]  Direction 

The  [ill]  direction  is  completely  non- degenerate  in  k  -space. 
That  is,  the  shadows  of  each  of  the  six  ellipsoids  do  not  overlap  for 
reasonably  low  energies.  However,  the  ellipsoids  are  degenerate  in 
energy,  and  so  the  general  formula,  equation  (55),  reduces  to: 

P(E)  -  [2ef(E)/4i  (2tt)3]  x 

a  b 

6  x  ( /(/d(E  -  E,  )dk,)dkt)  (54) 

-a  -b 


a  -  [2  m  E/ti2]1 

b 

b«  [  (2  m* /-II2)  (E  -  -fi2k  z2/2m,  J1 

(55) 

m#  »  (2m  +  m  )/3 

X  b 

k0  *  [4  (3 •  7 eV )  m*/5h2]l/2 

(56) 

E,  »  [-ft 2( k,  +  k0)2/2m*  +  il2k2/2m,] 

(57) 

Valence  Band  Calculation 

The  valence  band  calculation  for  P(E)  proceeds  from  a 
determination  of  the  transverse  energy.  The  calculation  is  simpler  from 
the  standpoint  that  the  valence  bands  are  maximal  in  energy  at  te’O  (i.e. 
hole  energy  is  minimal  at  k*0).  Thus,  the  large  transverse  momenta, 
associated  with  the  conduction  bands  of  indirect  gap  semiconductors  such 
as  silicon,  are  absent  here. 

We  start  with  the  general  form  of  the  constant  energy  surfaces 
for  the  valence  band.  This  equation  is  [8]: 

E  »  Ev  -  (-h2/2m)  x  [Ak2  +  [B2k4  ♦  C2(kx2kv2  +  kx2 kz2  + 
ky2kx2)]l/2]  (58) 

with  A  *  4.0  ;  B  —1.1  ;  C  —4.1  (59a) 

This  yields  three  effective  masses  for  the  valence 
band: 


ml  -  m/(A-B)  ■  m/2.9  ;  m^  “  m/(A+B)  ■  m/5.1 
and  m^  -  m/A  *  m/4  (59b) 

Taking  the  hole  energy  to  be  positive  downward  [47],  relabelled  E',  we 
get: 


P(S') 


3 

I  [2ef(E')/fi(2*)3]  x  (2w) 


l  =  i 


i 


z 
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•'mi 

y  B(E 1  -  E,  +  .V  ♦  Eg)  kt  dkt 

(60) 

0 

with 

kmi=  [2m,  E’/*2]1  72 

(61a) 

and 

Et  -  «2k,^2m, 

(.61  b) 

or,  rewri 

ting  in  terms  of  E'  ■  -  B  -  E,,  gives,  for  E  i 
Pvb  (E)  -  [2ef(E)/ft  (2*)3]  x  2*  x 

-V 

I  [  /  D(-E+-At?2m,  ♦  X)  k.  dkt] 

j=1  /  t  , 

(62) 

with 

a  -  [2m,(  -  E  -  Eg  i/H2]1  72 

(63) 

Predicted  Energy  Distributions  from  Si  [lOO] 

The  equations  developed  in  the  previous  sections  allow 
calculation  of  the  expected  FEEDs  from  Si.  In  particular,  predictions 
for  the  [too]  direction  appear  in  Figures  8  through  10.  The  effect  of 
raising  the  electron  temperature  is  demonstrated  by  the  broadening  of 
the  conduction  band  peak,  demonstrated  in  Figure  8.  Alternatively,  the 
effect  of  temperature  broadening  on  the  valence  band  emission  appears  in 
Figure  9.  By  raising  the  Ferai  level  above  the  conduction  band  edge, 
the  conduction  band  distribution  approaches  that  predicted  for  a  free 
electron  metal,  as  shown  in  Figure  10. 


In  general,  for  a  conduction  band  only  a  few  tenths  of  an  eV 
below  the  Fermi  level,  the  300K  distribution  consists  of  a  sharp  spike. 


.025 


Figure  8.  Effect  of  Temperature  on  (100)  Distributi 
Degenerately  N-Type  Si. 


0.4  V/A 
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Figure  9.  Effect  of  Temperature  on  (100)  Distribution  from 
Degenerately  P-Type  Si. 


28 


evidenced  in  Figure  10.  For  the  normalized  distributions  plotted,  the 
valence  band  will  be  of  low  intensity.  Of  course,  for  emission  direc¬ 
tions  other  than  the  <100>  set,  conduction  band  intensity  will  be 
greatly  reduced.  This  stems  from  the  non-zero  transverse  energy  terms 
within  the  Stratton  integrals. 

Conversely,  the  valence  band  distribution  does  not  change 
shape  if  the  Fermi  level  remains  above  it  in  energy.  In  such  a  case, 
the  valence  band  contribution  only  changes  amplitude.  However,  if  the 
Fermi  level  is  positioned  within  the  valence  band,  the  distribution  from 
this  band  narrows,  as  demonstrated  in  Figure  10,  where  both  cases  ap¬ 


pear 


CHAPTER  III 


APPARATUS  AND  EXPERIMENTAL  PROCEDURES 

Mechanical  Components 
Vacuum  System 

An  ultra-high  vacuum  (UHV)  system  is  employed  in  the  FE  study. 
The  UHV  system  is  comprised  of  a  stainless  steel  (ss.)  chamber,  with 
components  connected  to  its  Cu-sealed  flanges.  These  components  are 
viewports,  various  ceramic  feedthroughs,  valves,  a  "cold  finger”  sample 
holder,  and  a  nude  Bayard-Alpert  ionization  gauge.  In  addition,  a 
sealed  palladium  tube,  silver  brazed  to  a  miniConflat  flange,  is 
attached  to  the  system.  When  heated  with  an  alcohol  lamp,  the  palladium 
tube  diffuses  hydrogen  into  the  chamber. 

The  chamber  is  evacuated  by  various  standard  pumps.  From 
atmospheric  pressure,  the  system  is  initially  roughed  to  10-^  torr 
by  an  Edwards  Model  SD-1  00  rough  pump,  via  a  Varian  bakeable  mini¬ 
valve,  mounted  below  the  chamber.  The  rough  pump  is  then  switched  to 
back  an  Edwards  Model  E0-2  oil  vapor  diffusion  pump  (d.p.).  The  d.p. 
utilizes  Dow  Corning  DC-705  pump  fluid.  A  liquid  nitrogen  cold  trap, 
stationed  above  the  d.p.,  reduces  backstreaming  of  oil  vapor,  and  is 
connected  to  the  chamber  with  a  Granville-Phillips  1.5-inch  main  valve. 


Figure  11.  Schematic  Layout  of  Ultrahigh  Vacuum  System. 

Here,  i.g.  refers  to  the  3ayard-Alpert  ionization 
and  d.p.  refers  to  the  oil  vapor  diffusion  pump. 
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Also ,  attached  to  the  rear  of  the  chamber  is  a  titanium  filament 
sublimation  pump,  or  "getter,"  which  may  be  cryogenically  cooled.  A 
schematic  of  the  chamber  and  pumping  system  appears  in  Figure  11. 

After  the  chamber  pressure  drops  below  5  x  10“^  torr,  the 
system  is  baked  to  280°c  for  12  hours  under  a  fiberglass  oven. 

During  this  bakeout,  the  getter  is  "flashed"  by  passing  45  Amps  through 
the  filament  for  3  minutes.  This  effectively  degasses  the  filament  and 
getter  walls.  Periodically  during  baking,  the  getter  is  similarly 
flashed.  Finally,  the  heat  is  switched  off  and  the  system  is  cooled 
slowly  under  the  oven  to  reduce  thermal  strain  on  the  ceramic 
feedthroughs.  The  pumps  bring  the  ultimate  pressure  below  4  x 
10  11  torr  after  one  bakeout,  as  measured  with  the  integral  nude 
Bayard-Alpert  ionization  gauge. 

Energy  Analyzer 

A  high- resolution  van-Oostrum  type  electron  energy  analyzer, 
chosen  for  ease  of  construction  and  operation,  and  reproducibility  of 
results,  was  built  by  the  author.  Specifications,  published  by  Swanson 
and  Crouser  [53]*  were  modified  slightly  for  mounting  and  viewing  ease 
in  the  present  ’JHV  chamber.  An  Au-coated  Cu  cup  comprises  the  final 
collector;  the  rest  of  the  analyzer  is  constructed  of  ty  ^  304  ss., 
except  the  viewing  screen  composed  of  phosphor  coated  Ni  foil.  Present 
sizes  of  the  apertures  are  as  follows: 


(a)  Screen  probe  hole:  0.069  inch 


(b)  Lens:  0.139  inch 

(c)  Faraday  cup:  0.273  inch 
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The  analyzer  is  shown  in  cross  section  in  Figure  12. 

Though  a  lOmeV  resolution  was  claimed  by  design,  the  measured 
resolution  of  the  analyzer  was  only  slightly  better  than  HOmeV,  as 
judged  by  the  Young-Kuyatt  criterion  [59].  This  criterion  defines  the 
resolution  as  the  energy  difference  between  the  \Q%  and  30%  heights  of 
the  leading  edge  of  a  0°K  free  electron  emitter.  The  resolution  was 
not  improved  noticeably  by  the  inclusion  of  magnetic  shielding  within 
the  UHV  chamber.  However,  the  resolution  was  significantly  decreased  by 
the  presence  of  any  external  magnetic  field  sources  in  the  vicinity  of 
the  chamber;  for  this  reason,  no  ion  pump  was  installed  on  the  chamber. 

A  typical  tungsten  calibration  run  appears  in  Figure  13* 

One  of  the  changes  made  in  the  design  of  the  analyzer  was  to 
allow  the  collector  Faraday  cup  to  be  removed.  This  allows  optical 
alignment  of  the  analyzer  with  the  tip  of  the  sample.  The  Ni  foil 
screen  was  modified  from  a  flat  plate  into  a  semi- spherical  truncated 
cone,  to  reduce  distortion  of  the  electron  image. 


Tungsten  Calibration  Diatribution. 
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Total  Current-Voltage-Light  Measurements 

Although  the  energy  analysis  is  computer  interfaced,  total 
current  measurements  as  a  function  of  light  or  anode  voltage  are 
recorded  manually.  In  this  case,  light  is  supplied  to  the  right  side 
window,  normal  to  the  axis  of  the  sample,  by  a  Bausch  and  Lomb  High 
Intensity  Monochromator.  The  monochromator  is  comprised  of  a  model  BL33- 
86-39-01  Tungsten  (Quartz-Iodine)  Light  Source,  coupled  to  a  model  3L33- 
86-77  digital  readout  IR  Grating.  Slits  on  the  monochromator  are 
adjusted  for  a  spectral  width  of  19nm.  Two  filters  are  employed  to 
eliminate  the  second  order  from  the  grating--a  #2403  Pyrex  filter,  with 
cutoff  below  605nm,  and  a  #7-56  Corning  filter,  with  a  900nm  cutoff. 

With  these  filters,  the  operational  regime  of  the  monochromator  is  605nm 
to  1 800 nm. 

Anode  power  is  supplied  by  a  Computer  Power  Systems  model  CPS- 
100  0-30kV  regulated  power  supply.  The  output  from  this  supply  is 
applied  both  to  the  anode,  via  the  25kV  Ceramaseal  feedthrough,  and  a 
20QM:2Q0kohm  Victoreen  resistor  voltage  divider.  The  divided  anode 
voltage  is  measured  with  a  4.5  digit  Data  Precision  MM  model 
248. 

Before  taking  1-7  data,  the  surface  of  the  sample  is  observed. 


This  is  undertaken  by  applying  about  IkV  to  the  phosphor  screen  of  the 
analyzer.  Screen  voltage  is  provided  by  a  Fluke  model  408B  High  Voltage 
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maintained  at  about  100-200  volts.  At  this  potential,  negligible  light 
intensity  is  generated  by  the  phosphor,  eliminating  any  photo¬ 
stimulation  feedback  to  the  sample. 

Measurement  of  the  surface  potential  at  constant  current  is 
obtained  in  the  following  way.  As  shown  in  Figure  14,  the  collector  cup 
is  maintained  at  virtual  ground  using  the  Current- to-Voltage  Converter 
(CVC),  described  in  the  Electronics  section  following,  to  measure  the 
probe  hole  current.  Output  from  the  CVC  is  monitored  on  a  Data 
Precision  DVM.  The  anode  potential  is  fixed  in  positive  value.  The 
negative  tip  potential  is  varied  until  a  current  increase  is  observed. 
This  value  is  deemed  the  "cutoff  potential."  The  tip  potential  is 
measured  simultaneously  using  another  Data  Precision  DVM.  The  cutoff 
potential,  anode  voltage,  and  probe  hole  current  readings  are  manually 
recorded.  Following  these  measurements,  the  total  current  is  measured 
by  applying  the  sum  of  the  previous  anode  and  offet  voltages  to  the 
anode  ring,  and  connecting  an  ammeter  to  the  sample.  Again,  the  sample 
is  thus  at  "virtual  ground"  potential. 


Surface  Potential  Measurements 


Electronics 


Computer  Interfacing 

Hardware 

The  energy  analyzer  is  completely  run  by  a  microcomputer. 

This  provides  rapid  data  acquisition,  uniform  and  precise  measurement, 
signal  averaging,  numerical  differentiation  and  smoothing,  and  data 
storage.  Further,  both  CRT  and  plotter  outputs  are  available  for  rapid 
viewing  or  hardcopy  output  of  data. 

Outlined  in  Figure  15  is  a  flow  chart  of  the  hardware 
interfacing.  Beginning  with  the  Apple  II  microcomputer  (uP),  decoding 
of  specific  address  codes  sets  an  initial  potential  on  the  retarder 
supply.  The  current- to- voltage  converter  (CVC)  measures  the  resulting 
current  on  the  retarder  cup.  This  signal  is  passed  through  a 
Butterworth  filter  and  subsequently  measured  by  the  ADC,  which  converts 
it  to  a  12-bit  binary  number.  In  two  operations,  the  uP  reads  the 
number,  and  adds  it  to  previously  accumulated  data  for  the  particular 
energy  level. 

Various  addresses  are  reserved  for  special  purposes.  These 
address  codes  are  decoded  outside  the  computer  and  utilized  to  switch 
counters,  control  lines,  or,  in  conjunction  with  the  data  bus,  to  put 
digital  values  onto  the  plotter  digital- to-analog  converters  (DACs). 
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TABLE  1 

Hardware  Addresses  and  Data  Assignments. 


Address 

Assignment 

Decimal 

Hexadecimal 

28680 

$7007 

V-Min,  Byte  *4,  MSB 

28681 

7008 

V-Min  ♦  2. 54mV ,  Byte  #4,  MSB 

28632 

7009 

V-Min  ♦  5 . 08mV ,  Byte  #4,  MSB 

32768 

8000 

R,  #Readings/2. 54mV  step 

32769 

8001 

L,  #Loops  of  total  voltage  range 

32770 

8002 

YA 

32771 

8003 

YB 

32772 

8004 

GA 

32773 

8005 

GB 

32774 

8006 

VI  ■  Tip  voltage,  IB 

32775 

8007 

VI,  HB 

32776 

8008 

Anode  voltage,  IB 

32777 

8009 

Anode  voltage,  HB 

32778 

800A 

Tip  type 

32779 

800B 

Month 

32780 

800C 

Day 

32781 

800D 

Year 

32782 

800 E 

32783 

800F 

Lens  voltage/2 

32784 

8010 

Temperature  in  Kelvin/4 

32785 

8011 

Offset  voltage,  IB 

32786 

8012 

Offset  voltage,  HB 

32787 

8013 

32738 

8014 

32789 

8015 

V-min,  IB 

32790 

8016 

V-min,  MB 

32791 

8017 

V-min,  HB 

32792 

8018 

V-min  *  2. 54 mV,  IB 

32793 

8019 

V-min  +  2.54mV,  MB 

32793 

801 A 

V-min  +  2.54mV,  HB 

32794 

e 

8013 

V-min  ♦  5.08mV,  IB 

Notes:  For  data  storage,  S7007  to  S79FF  and  S8015  to  $BFFF: 

IB  *  Lowest  Byte  MB  *  Next  Highest  Byte  KB  *  Next  Highest  Byte 
MSB  *  Most  Significant  Byte 


For  parameters,  at  the  remaining  locations  in  this  table: 
12  *  Low  3yte  H3  ■  High  Byte 
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TABLE  1  (continued) 

Hardware  Addresses  and  Bata  Assignments. 


Address 

Assignment 

Becimal 

Hexadecimal 

49410 

SC102 

Increment  Retarder  DAC  by  2.54meV. 

49412 

Cl  04 

Reset  Retarder  to  lOeV 

49416 

C108 

Data  13 

49417 

C109 

Data  High  Nybble 

49418 

C10A 

Reset  and  Start  ADC 

49419 

Cl  OB 

Sample  and  Hold:  Hold 

49420 

C10C 

Sample  and  Hold:  Sample 

49421 

C10D 

Data  Bus  Diagnostic,  ■  255 

49422 

C10E 

Data  Bus  Diagnostic,  *  0 

49424 

C110 

North  Plotter  Motion,  2.54eV  increment  t 

49425 

cm 

Northeast  Motion  * 

49426 

Cl  1 2 

East  Motion  -*■ 

49427 

Cl  1 3 

Southeast  Motion  > 

49428 

Cl  14 

South  Motion  ♦ 

49429 

Cl  1 5 

Southwest  Motion  ^ 

49430 

Cl  1 6 

West  Motion 

49431 

Cl  1 7 

Northwest  Motion  H 

49432 

Cl  1 8 

Pen  Up— Not  Connected 

49433 

Cl  1 9 

Pen  Down— Not  Connected 

49434 

Cl  1 A 

Load  X  HB 

49435 

Cl  IB 

Load  X  Low  Nybble 

49436 

Cl  1 C 

Load  Y  KB 

49437 

Cl  ID 

Load  Y  Low  Nybble 

49438 

cue 

Clear  X 

49439 

Cl  IF 

Clear  Y 
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The  address  codes  used  in  the  present  interfacing  are  presented  in  Table 
1 . 

The  sections  of  the  block  diagram  are  shown  in  Figures  16 
through  21.  Each  section  is  described  briefly  below. 

1.  Buffers.  This  wire-wrapped  board,  shown  in  Figure  16, 
slips  into  the  Apple' s  interface  slot  #1  .  It  protects  the  Apple  from 
outside  transients,  and  is  powered  from  the  Apple's  own  power  supply. 
This  board  contains  LS367  chips  as  unidirectional  address  bus  buffers, 
and  a  National  Semiconductor  #8304  octal  data  bus  transceiver. 

2.  Decoding  and  Plotter  Control.  This  wire-wrapped  board, 
shown  in  Figures  17  and  18,  preliminarily  reduces  the  designated  16-bit 
addresses  to  managable  4-bit  addresses.  The  154  chip,  shown  in  Figure 

17,  accomplishes  this  decoding.  The  addresses  go  to  specific 
peripherals  called  by  the  uP  program.  In  addition,  certain  addresses 
are  completely  decoded  (using  chips  numbered  154  and  LS42)  to  single 
commands,  connected  to  cascaded  4-bit  counters  (LS193  chips).  These 
counters  drive  two  12-bit  DACs  (Analog  Devices  AD  DAC80-CBIV)  for  the 
"X"  and  “T"  directions  of  the  plotter.  The  data  bus  allows  any  preset 
number  to  be  placed  on  these  DACs.  This  circuitry  is  depicted  in  Figure 

18. 

3.  Retarder  Controller.  Shown  in  Figure  19,  this  section 
consists  of  a  transformer  isolated  power  supply  and  7.5kV  optoisolators, 
driving  an  isolated  12-bit  DAC  (Analog  Devices  AD  DAC80-C3IV).  This  DAC 
is  connected  to  provide  precise  2.54mV  potential  3teps  to  the  tip 
through  a  unity  gain  741  Op  Amp  buffer.  Only  two  controls  are  utilized- 
-a  preset  to  -10V,  and  an  increment  command.  The  next  increment  after  a 
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gure  15.  Block  Diagram  of  Electronics 


Figure  16.  Buffer  Board  for  Apple  II  Interfacing. 
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Figure  20.  Current- to-Voltage  Converter  (CVC)  and  Butteruorth 
Filter. 
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Figure  21.  Analog- to-Digital  Converter  (ADC). 
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preset  command  sets  the  DAC  to  zero  volts. 

4.  Current- to-Vol tag e  Converter  (CVC).  The  CVC  consists  of  an 
ultra- high  impedance  (10^  ohm)  precision  FET  Op  Amp  (Analog  Devices 
AD545KD)  connected  in  the  current- to- voltage  converter  mode.  Mounted  on 
a  solid  teflon  sustrate  to  minimize  leakage  currents,  this  CVC  features 
interchangeable  feedback  resistors  R  ,  giving  output  voltages  of  -I  x 
R  .  The  output  from  the  AD545  is  connected  to  a  Butterworth  low  pass 
filter.  The  CVC  is  illustrated  in  Figure  20. 

5.  Analog- to-Digital  Converter  (ADC).  The  ADC,  shown  in 
Figure  21,  is  a  12-bit  chip  (Analog  Devices  AD574KD),  with  a  sample-and- 
hold  (SdH)  on  its  input.  The  S4H  (Analog  Devices  AD582KD)  maintains  the 
input  signal  at  a  constant  value  during  the  "read"  cycle  of  the  AD574. 

Ia  front  of  the  SJtH  are  connected  three  FET  input  Op  Amps  in  an 
instrumentation  amplifier  scheme,  to  eliminate  the  common  mode  present 
in  the  ground  loop  from  the  UHV  chamber  to  the  ADC. 

6.  Offset.  Offset  voltage  is  provided  by  a  Hewlett-Packard 
model  HP6110A  regulated  DC  power  supply.  The  offset  compensates  for 
potential  differences  between  the  sample  surface  and  ground.  The  offset 
need  be  applied  only  if  the  energy  distribution  cannot  be  measured  with 
the  lOeV  range  attainable  with  the  bare  DAC.  Otherwise  this  point  is 
grounded . 


Software 

The  controlling  program  for  data  acquisition  has  two  distinct 


parts.  The  first,  programmed  in  Basic,  allows  the  user  to  input 
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essential  experimental  information  to  the  computer.  This  information  is 
requested  by  the  program  input  statements,  and  stored  along  with  the 
actual  acquired  data  for  future  reference.  The  date,  sample  type,  lens 
and  anode  potentials,  offset  voltage,  and  sample  temperature  are  polled, 
as  well  as  the  number  of  readings  (r)  and  loops  (l»)  desired  and  starting 
and  ending  voltages  for  each  loop.  These  last  two  parameters  are  called 
YA ,  YB ,  GA ,  and  G8  in  Table  1 . 

The  second  part  of  the  control  program  is  in  6502  machine 
language.  Though  less  interactive  with  the  user,  the  machine  language 
performs  1000  times  faster,  and  thus  provides  the  necessary  data 
acquisition  speed.  Multiple  precision  addition  can  also  be  performed 
easily  and  rapidly  in  machine  language.  The  machine  program’ s  quadruple 
precision  ability  allows  for  2 bits  of  information  to  be  stored 
for  each  2. 54meV  energy  increment,  corresponding  to  22®  full  scale 
readings.  This  capability  remains  unused;  for  most  applications, 
setting  the  experiment  through  L*50  loops  and  R“50  readings  suffices. 
However,  for  increased  data  smoothing  or  large  feedback  resistances, 
which  require  longer  settling  times,  up  to  R*255  and/or  L«255  are 
possible. 

The  actual  programs  are  contained  on  cassette  tapes  in  the 
laboratory,  and  in  an  appendix  to  the  thesis.  One  merely  performs  a 
]lOAD  from  the  Basic  program  into  the  Apple,  then  flips  the  tape,  types 
]CALL  -151  to  put  the  computer  into  a  machine  language  mode,  then  types 
•6B00.6CFFR  for  reading  the  tape  into  the  proper  memory  locations.  The 
user  holds  the  Control  button  down  and  hits  C  to  return  to  the  Basic 
program,  then  types  ]RUB  to  use  the  program. 
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Although  most  elements  of  the  program  are  readily  prompted  by 
the  3asic  program,  the  plotting  software  i3  slightly  complicated.  The 
complications  arise  from  the  fact  that  the  accumulated  data  appears  as 
l(V),  and  what  is  desired  is  dl/dV.  Though  this  is  a  simple  matter  to 
take  dl/dV  *  [l(i+l)  -  l(i-l)]/2,  any  fluctuations  appearing  in  the 
original  data  cause  large  oscillations  in  the  calculated  derivative.  To 
smooth  the  oscillations,  one  can  change  the  definition  of  dl/dV  to 
dl/dV  -  [i  (i+m)  -  I  (i-m)  ]/(2m) 

which  reduces  the  oscillations  by  the  factor  m.  The  author  has  found 
that  another  manipulation  will  reduce  the  oscillations  even  further. 

This  last  method  is  to  average  the  values  at  the  points  (i+m),  and  (i- 
m) .  Functionally,  this  amounts  to  defining  l(i+m)  to  be 

n 

l(i+m)  -  [  I  l(i+nt+ j)]/[2n+1  ] 
i=-n 

This  implies  that  dl/dV  should  be 

n 

dl/dV  -  (  £  [l(i+m+ j)  -  l(i-m+ j)])  /  [(2n+l)s2m] 
j  =  -n 

By  appropriately  choosing  "m,"  called  the  "points  to  be  averaged,"  and 
"n,"  the  "mesh,"  the  data  may  be  smoothed  very  reasonably  without 
eliminating  any  inherent  characteristics  of  the  energy  distributions. 


Sample  Preparation 


Cutting  Blanks  from  Vafers 

Two  methods  have  been  successfully  employed  for  creation  of 
sample  blanks  in  the  "bridge  structure"  of  D'Asaro  [14].  The  first, 
useful  for  thin  wafers,  of  about  10  mils  thickness,  utilizes 
photolithography  and  produces  specimens  of  precise  dimensions.  The 
second  method,  useful  for  wafers  of  about  40  mils  thickness,  relies  upon 
hand  cutting  on  an  abrasive  disk,  glass  cutting  saw.  These  two  methods 
are  described  below. 

Before  performing  either  method,  the  desired  crystal 
orientation  is  first  determined.  This  may  be  done  by  observing  the 
backscattered  Laue  X-ray  pattern,  or  by  utilizing  the  cleavage 
characteristics  of  the  sample.  It  is  known,  for  example,  that  Si 
cleaves  exclusively  along  the  (ill)  plane.  If  one  is  using  <100> 
oriented  wafers,  the  <110>  direction  lies  normal  to  the  edge  of  the 
cleaved  wafer,  and  is  coplanar  with  the  wafer  surface.  This  is 
illustrated  in  Figure  22.  This  procedure  eliminates  the  need  to  X-ray 
the  samples  if  <100>  or  <110>  orientations  are  desired. 
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Photolithography 

Photolithography  is  a  technique  used  extensively  in  solid 
state  device  manufacturing  for  precisely  aligning  and  defining  circuit 
elements,  using  an  HP-acid  resistant  photoresist.  The  author  has 
adapted  this  technique  for  the  manufacture  of  thin,  precisely  defined  Si 
sample  blanks  for  use  in  FEM  work. 

Systematically,  the  method  consists  of  first  coating  both 
sides  of  each  wafer  with  Shipley  AZ-1350  positive  resist.  The  wafers 
are  baked  overnight  at  90°c  in  a  dry  nitrogen  ambient  for  good 
adhesion  to  the  surface.  Next,  the  wafers  are  clamped  in  a  jig  with  two 
mirror- image  contact  masks,  shown  in  Figure  23.  Both  sides  of  each 
wafer  are  exposed  through  the  masks  to  a  high  intensity  UV  source  for 
one  to  two  minutes,  depending  upon  photoresist  thickness.  Thereon,  the 
exposed  wafers  are  developed  in  corresponding  Shipley  developer,  rinsed 
in  deionized  water  (D.I.),  and  spun  dry.  The  wafers  are  rebaked 
overnight  as  before. 

Under  a  protective  hood,  one  to  four  wafers  are  then  held 
vertically  in  a  slotted  teflon  basket.  They  are  immersed  into  a  plastic 
beaker  of  48?  HP,  just  large  enough  to  contain  and  cover  the  wafers. 
Slowly,  several  drops  of  nitric  acid  are  mixed  into  the  beaker.  As 
these  acids  react  with  the  uncoated  Si  surface,  gas  bubbles  are  evolved 
and  the  surface  of  the  Si  slowly  dissociates.  As  the  surface  etches, 
however,  it  heats,  speeding  the  reaction.  Thus,  it  is  important  to  keep 
the  wafers  immersed  in  the  liquid  etchant  bath  else  the  resist  will  peel 
from  the  surface.  In  spite  of  the  best  precautions,  the  resist  may 


Figure  23.  Photolithographic  Mask  Designs. 

(A)  Photograph  of  Hubylith  master. 

(3)  Contact  prints  of  front  and  back  masks 
reduced  20X  from  master. 


begin  to  peel  from  the  wafers.  If  this  occurs,  the  wafers  must  be 
removed,  cleaned,  and  recoated  with  resist. 

When  etching  is  complete,  samples  float  to  the  surface  of  the 
acid  bath,  and  are  retrieved  with  a  pair  of  tweezers.  Subsequently, 
they  are  rinsed  in  D.I.  water,  then  in  acetone  to  remove  the  Shipley 
resist,  and  finally  in  D.I.  water. 

These  samples  are  now  sandblasted  gently  to  round  the  tip  end 
After  this  treatment,  the  samples  are  etched  electrolytically,  as 
described  in  the  next  section. 

Abrasive  Cutting 

This  technique  provides  roughly  equivalent,  rapidly  generated 
samples.  It  requires  wafers  of  standard  4-inch  diameter  thickness, 
about  40  mils.  Typically,  one  such  4-inch  wafer  is  scribed  and  broken 
into  quarters.  Then  these  are  cemented  onto  a  3/1  6- inch  thick  glass 
substrate  with  melted  Cenco  red  wax.  The  glass  and  four  stacked 
quarters  are  cut  by  hand  on  a  standard  glass- cutting  abrasive  wheel  saw 
The  glass  substrates  and  Si  blanks  are  heated  and  the  wax  blotted  off, 
then  the  blanks  are  cleaned  in  acetone.  Finally,  the  tip  ends  of  the 
blanks  are  sandblasted  round.  The  rounding  is  necessary  because  the 
etching  action  is  isotropic--elliptical  cross  sections  are  etched  into 
"butter  knives"  with  one  dimension  sharp  and  the  other  blunt. 
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Electrolytic  Etching 

Electrolytic  etching,  rather  than  purely  chemical  etching, 
must  be  used  for  creation  of  successful,  utilizable  tip3.  Best  results 
are  obtained  using  the  techniques  described  by  Busch  and  Fischer  [ll], 
who  suggest  optimum  ratios  of  HFsHITO^  for  various  dopings  of  Si.  In 
general,  the  ratio  is  5s 5  for  resistivities  greater  than  10  ohm-cm,  and 
drops  to  Is  10  for  resistivities  of  0.1  ohm-cm  or  less.  Equally 
important  is  the  use  of  a  variable  AC  voltage  supply,  with  a  platinum 
wire  immersed  in  the  acid  bath  and  another  electrode  attached  to  the  Si 
sample.  The  voltage  is  regulated  until  small  bubbles  appear  at  both 
sample  and  Pt  wire  electrodes.  This  voltage  may  vary  from  a  few  tenths 
of  a  volt  for  high  conductivity  samples  to  60  or  more  volts  for  nearly 
intrinsic  samples.  Lastly,  the  sample  is  rinsed  in  D. I.  water,  dried 
with  compressed  nitrogen,  and  stored  in  a  closed  container. 

The  electrode- to-Si  contact  must  be  excellent  to  assure 
consistent  results.  To  this  end,  the  author  always  strips  the  oxide 
from  the  "legs"  of  the  sample  with  etchant,  rinses  the  sample,  then 
attaches  a  strong,  flat  clip.  Loss  of  etching  activity  is  almost  always 
attributable  to  poor  electrode- to- sample  contact. 
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!'  Prevacuum  Processing 

i 

i  Just  prior  to  installation  in  the  UHV  system,  a  previously 

etched  tip  is  removed  from  the  closed  container  and  given  an  "RCA 
clean."  This  consists  of  first  a  4-minute  dip  in  boiling  NH^OH: 

(4:10:1),  followed  by  a  rinse  in  D. I.  water.  Secondly,  the 
sample  is  dipped  into  boiling  HC1:H20:H202  (4:10:1),  and  again  rinsed 
in  D. I.  water.  Lastly,  the  sample  is  dipped  in  48%  HP,  rinsed  twice, 
and  dried.  After  insertion  into  the  cold  finger  tip  holder,  a 
resistance  measurement  between  the  tip  feedthroughs  indicates  the 
contact  resistance.  If  considered  too  high,  another  dip  in  HF  and  rinse 
usually  eliminates  the  problem. 

In-Vacuum  Processing 

i 

Once  the  sample  has  been  mounted  in  the  chamber,  the  system  is 
pumped  out.  The  procedure  for  evacuation  is  described  in  the  section 
entitled  "Vacuum  System."  After  low  10”^ ^  torr  has  been  obtained  after 
bakeout,  the  experimenter  is  able  to  clean  the  sample  with  the  following 
technique. 

With  pure  hydrogen  admitted  through  the  sealed  Pd  tube  to  a 
pressure  of  10”^  torr,  a  novel  method  is  employed  to  clean  the  Si. 
Following  the  work  of  Busch  and  Fischer  [ll],  who  described  a  faintly 
,  visible  DC  field  desorption  image  of  Si  in  H  ,  the  author  performs 
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an  AC  field  desorption- field  emission  cleaning.  This  method  allows  a 
better  view  of  the  surface  cleaning  than  the  DC  desorption  technique. 
Following  Cooper  and  Muller  [l 3 ] *  who  first  used  the  technique  to  clean 
tungsten  in  vacuum,  an  adjustable  high  voltage  AC  is  applied  to  the  tip 
while  an  adjustable  negative  DC  voltage  is  applied  to  the  ring  "anode" 
of  the  analyzer.  A  positive  voltage  of  3  kV  is  put  on  the  phosphor 
screen.  The  resulting  field  at  the  sample  surface  oscillates  between 
field  emission  and  field  desorption.  The  field  emission  gives  the 
experimenter  an  excellent  view  of  the  sample  surface  as  the  cleaning  is 
occurring.  The  reacts  with  the  silicon  dioxide  during  desorp¬ 
tion,  enhancing  the  removal  of  this  strongly  bonded  impurity,  and 
exposes  the  silicon  surface. 

While  the  AC  field  is  on,  the  i-8  pumped  out.  After  a 
short  time,  the  AC  is  cut  off,  and  the  system  is  again  baked.  This 
second  bakeout  leaves  the  system  in  low  10“^ 1  torr.  At  this  point, 
another  AC  field  desorption- field  emission  cleaning  is  performed,  this 
time  in  vacuum  without  hydrogen. 

If  the  "RCA  clean"  procedure,  outlined  in  the  section  on 
prevacuum  processing,  is  followed  rigorously,  the  author  has  found  that 
the  hydrogen  reaction  described  above  may  be  eliminated  altogether. 

This  allows  the  experimenter  to  maintain  a  sharp  tip  from  the  onset,  and 
obviates  the  extra  bakeout  required  with  the  hydrogen  cleaning. 

If  annealing  is  desired,  a  current  is  passed  through  the  tip 
"legs,"  and  the  sample  is  resistively  heated  and  annealed.  As  Si  has  a 
negative  coefficient  of  resistance,  ballast  protection  resistors  must  be 
used  for  high  resistivity  samples  [43].  For  low  resistivity  samples,  a 


transistor- regulated  cons taut- current  power  supply  suffices  to  heat  the 
sample,  automatically  adjusting  for  resistance  changes. 
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CHAPTEH  IV 

EXPERIMENTAL  RESULTS 

Results  of  the  experimental  measurement  of  the  field  emission 
energy  distribution  (PEED),  and  Fowler  Nordheim  current- voltage  (F-N) 
plots,  from  silicon  of  various  dopings  are  presented  in  this  section. 

It  was  observed  by  the  author  that  surface  conditions,  as  well  as  tem¬ 
perature  and  doping,  affected  the  results.  For  this  reason,  photo¬ 
graphs  of  the  various  surface  conditons  are  presented  in  Figure  24. 

Lightly  Doped  Silicon 

Energy  Distributions 

One  of  the  main  objectives  of  the  thesis  was  to  ascertain  ob¬ 
vious  differences  in  the  FEED  of  Si  under  PFE  conditions.  To  this  end, 
samples  were  cut  from  very  pure  p-type  Si  wafers.  These  samples  yielded 
results  typified  in  Figure  25 •  The  FEEDs  obtained  from  these  samples 
were  characteristically  wide  (FVHM  >  0.75eV)  and  Gaussian  in  shape  at 
both  300K  and  78K.  Further,  the  FEED  generally  contained  large,  non- 
reproducible  fluctuations  in  amplitude  at  points  of  the  distribution. 

The  only  obvious  difference  in  the  FEED,  between  light  and  dark 
conditions,  is  the  offset  potential.  There  are  no  distinguishing 
differences  in  the  shapes  of  the  distributions,  even  though  the  PFE 
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Figure  24.  Various  Surface  Conditions  Observed  for  Emitters. 

(a)  N-type,  0.001  ohm-cm,  field  desorbed, 

(110)  oriented,  side  view. 

(b)  Same  as  (a),  rear  view. 

(c)  P-type,  60-100  ohm-cm,  thermally  annealed, 
(110)  oriented,  side  view. 

(d)  Tungsten,  thermally  annealed,  (110)  oriented, 
side  view. 
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Figure  24  ( continued) . 

(e)  3-type,  0.001  ohm-cm,  thermally  annealed, 

(110)  oriented,  aide  view. 

(f)  Same  as  ( e) ,  rear  view. 

(g)  P-type,  60-100  ohm- cm,  field  desorbed, 

(110)  oriented,  rear  view. 

(h)  Same  as  (g),  thermally  annealed  at  T<700C, 

(110)  oriented,  rear  view,  showing  ring  formation. 
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figure  25.  Illwinatftd  vs.  Carle  FEEDs  from  Kigb  Resistivity 
P-Type  Si. 
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Pigure  26,  Photo-Field-Eaiasion  (PPE)  Fowltr-Nordheia  Plot 


TABLE  2 

Photo-Field  Emission  Surface  Potential  Measurements 


Anode-Surface 
Voltage  (Volts) 

— 

Base  Contact-Surface 
Voltage  (Volts) 

Total  Emitter  Current 
(Amps  x  10+1®) 

Dark 

Light 

Dark 

Light 

1650 

7.1 

6.7 

2.40 

2.41 

1680 

9.5 

8.7 

3.30 

3.40 

1710 

12.0 

11.0 

4.18 

4.35 

1740 

19.0 

17.0 

4.90 

5.20 

1770 

25.5 

23.7 

5.50 

5.95 

1800 

32.6 

31.5 

6.00 

6.40 

1830 

41.5 

39.9 

6.41 

6.85 

1860 

49.0 

47.6 

6.80 

7.23 

1890 

58.6 

56.9 

7.20 

7.70 

Scale 

Change 

(Amps  x 

10*9) 

0.75 

0.81 

1920 

68.0 

65.8 

•0.77 

0.84 

1950 

76.0 

73.8 

0.82 

0.87 

1980 

85.7 

85.7 

0.85 

0.92 

2010 

93.9 

92.4 

0.87 

0.94 

2040 

0.92 

0.97 

2070 

109.5 

107.4 

0.94 

1.00 

2100 

118.2 

116.3 

0.97 

1.04 

2130 

127.0 

125.7 

1.00 

1.06 

2160 

138.3 

135.3 

1.03 

1.08 

2190 

146.0 

144.0 

1.05 

1.12 

Notes:  Light  source:  B«L  High  Intensity  Monochromator 
Setting:  lOOOnm. 

Filter:  Pyrex  # 2403* 

Sample:  60-100  ohm-ca  p-type. 

Temperature:  JOCK. 

Surface:  Field  desorbed. 

Pressure:  2  i  10"  '  torr. 
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current  at  78K,  with  photo- stimulation,  was  usually  several  orders  of 
magnitude  larger  than  the  dark  current. 

Using  the  FEED  to  determine  the  offset  potential,  the  author 

measured  the  I-V  ^  relation  of  the  sample.  At  various  fixed 
offset  r 

offset  potentials,  the  PFE  current  was  always  slightly  higher  than  the 
dark  current  at  the  same  offset  potential.  Data  on  this  effect  is 
presented  in  Table  2. 

This  FEED  information  may  indicate  several  possibilities.  The 
wide  energy  distributions,  accompanied  by  fluctuations,  appear  to 
indicate  some  non-equilibrium  situation.  Indeed,  various  researchers 
have  argued  that  the  width  is  indicative  of  hot  electron  effects  in  the 
sample.  Distributions  up  to  4eV  in  width  were  reported  by  Kagan,  et  al. 
[32].  As  pointed  out  by  Stratton  [52],  this  implies  that  these  elec¬ 
trons  have  energy  larger  than  the  barrier  height  and  that  a  complete 
reformulation  of  field  emission  theory  is  then  necessary. 

Current-Voltage-Light  Characteristics 

Fowler-Nordheim  plots  of  data  from  lightly  doped  p-type  Si 
follow  the  same  pattern  regardless  of  surface  condition.  An  example  is 
shown  in  Figure  26.  However,  the  ratio  R  of  photostimulated  current  to 
dark  current  at  30CK  is  highest  for  the  field  desorbed  surface.  The 
effect  of  slight  annealing  at  T<700C  is  the  formation  of  a  ring  on  the 
previously  desorbed  surface,  and  a  consequent  reduction  in  R,  from  1.68 
to  1.42.  This  surface  is  shown  in  Figure  24g.  The  Current-Voltage 
plots  illustrating  this  effect  appear  in  Figure  27. 


field  oesorbed surface 

P-TYPE  SI  80-100  A-CM 
30  0'K 


CURRENT  iamperes 1 1 
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Also,  the  offset  potentials,  or  IR  drops,  indicated  non-ohmic 

1 2 

behavior,  with  resistances  of  up  to  10  ohms  at  300K,  as  indicated 
in  Figure  28. 

Additionally,  the  FFE  increase  varies  with  wavelength  and 
field.  This  behavior  is  demonstrated  in  Figures  29  and  30,  with  Log  of 
Current  vs.  Illumination  Wavelength  at  fixed  values  of  V  .  As  indi¬ 
cated,  the  photocurrent  is  substantially  larger  at  fixed  for  photon 
energies  larger  than  the  thermal  band  gap.  Nevertheless,  residual  en¬ 
hancement  does  exist  at  78K  for  photons  of  energy  down  to  1.03eV, 
corresponding  to  1200nm.  As  discussed  by  Borzyak,  et  al.  [lO],  this 
residual  enhancement  may  be  due  to  absorption  by  carriers  in  surface 
states,  or  perhaps  by  free  carrier  or  impurity  absorption. 

Heavily  N-Doped  Silicon 

Energy  Distributions 

Because  the  FEED  studies  on  lightly  doped  Si  showed  no  obvious 
differences  in  energy  structure,  the  author  decided  to  elicit  informa¬ 
tion  from  more  heavily  doped  samples.  It  was  hoped  that  the  increased 
doping  would  eliminate  the  large  offset  voltages  and  current  fluctua¬ 
tions  accompanying  the  previous  data.  Further,  it  was  expected  that  the 

19  3 

heavily  n-doped  Si  used  (0.001  ohm-cm  resistivity,  10  5  per  cm  Sb 
doped)  would  show  a  metallic  energy  distribution  and  provide  a  starting 
point  for  further  data  comparison. 

These  expectations  were  largely  unfulfilled.  Although  current 
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fluctuations  were  eliminated  for  desorbed  or  annealed  samples,  the  zero 
offset  distributions  expected  for  metallic  emitters  were  never  observed. 
Attempts  to  connect  the  non-zero  offset  voltage  to  an  IR  drop  were  com¬ 
pletely  unsuccessful.  Indeed,  the  "onset"  potential,  defined  by  the 
author  to  be  the  offset  voltage  yielding  10-1^  Amp  probe  hole  current, 
actually  decreased  with  increasing  total  tip  current.  This  behavior  is 
shown  with  the  current- voltage  data,  in  Figure  31. 

Field  Desorbed  Surface 

The  field  desorbed  surface,  shown  in  Figure  24b,  yielded 
reproducible  FEED  results.  These  results  appear  in  Figures  32  and  33* 
The  energy  distributions  may  be  characterized  as  follows  [28]: 

(a)  300K:  A  single  peak,  which  shifts  downward  in  energy  and 
broadens  with  increasing  field.  The  peak  is  initially  about  0.4eV  below 
Fermi  level. 

(b)  78K:  The  single  peak  broadens  and  shifts  as  in  (a),  but  at 
some  particular  field  there  arises  a  subsidiary  peak  from  below  the 
Fermi  level.  The  subsidiary  peak  increases  in  magnitude  with  respect  to 
the  larger,  low  energy  peak. 

Thermally  Annealed  Surface 

Thermally  annealed  surfaces  show  little  general  qualitative 
differences  to  the  results  from  desorbed  surfaces  at  300K  in  that  both 
types  show  an  energy  shift.  FEEDS  performed  on  the  surface  shown  in 
Figure  24e  yielded  a  single  peak,  which  shifted  down  in  energy  with 
increasing  field.  The  FEED  measurements  made  from  the  annealed  surface, 
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Figure  34.  300K  FEED  of  Heavily  N-Doped,  Annealed  Si 
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35«  Fovler-Uordheim  Plot  for  Heavily  N-Doped  Si 
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shown  in  Figure  34,  evidence  this  observation.  However,  the  general 
shape  of  the  energy  distribution  from  the  annealed  surface  appears 
slightly  different  from  that  of  the  desorbed  surface.  In  particular, 
the  annealed  surface  FEED  shows  a  longer  "tail"  on  the  lower  energy  edge 
of  the  distribution.  Such  a  change  in  the  distribution  shape  is  not 
predicted  for  valence  band  emission  in  the  Stratton  theory. 

Current-Voltage  Characteristics 

The  current-voltage  behavior  of  the  heavily  n-doped  Si  appears 
to  follow  the  F-N  equation  fairly  well.  A  F-N  plot  of  data  from  the 
field  desorbed  surface  appears  in  Figure  35.  No  severe  deviations  from 
F-N  behavior  are  observed,  in  contrast  to  the  data  obtained  from  high 
resistivity  p-type  material. 


Heavily  P-Doped  Silicon 

Only  one  sample  of  ( 1 1 0 )  oriented  P+  Si  was  studied.  This  >■ 

sample  had  a  resistivity  of  0.015  ohm-cm,  and  was  produced  from  a 
Monsanto- supplied  wafer.  The  surface  was  field  desorbed.  FEEDs  from 
this  sample,  which  appear  in  Figure  36,  show  a  single  peak,  located  1 

below  the  Ferai  level,  similar  to  the  data  obtained  from  heavily  n-  \ 

doped  Si.  Unlike  the  behavior  of  the  n-type  samples,  however,  the 
"onset"  potential  increased  with  increasing  field  as  demonstrated  in 
Figure  37.  Thus,  although  the  peak  appears  below  E  ,  it  also  may  be 
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Figure  36.  300 K  FEEDs  of  Heavily  P-Doped  Si  at  Various  Fields. 
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true  that  a  potential  drop  exists  within  the  sample.  The  F-N  plot 
associated  with  this  data  appears  in  Figure  38. 


Discussion  and  Comparison  of  Results 

Previous  Work  on  Si 

FEED  studies  have  been  performed  on  Si  and  other  semicon¬ 
ductors  previously.  For  Si,  Russell  and  Litov  [45]  reported  a  single 
peak  which,  upon  "fracture,"  two  peaks,  with  maxima  separated  by  the 
thermal  gap  energy  of  about  1.2eV.  There  appeared  to  those  authors  a 
simple  identification  of  the  lower  energy  peak  with  the  valence  band, 
and  of  the  upper  peak  with  the  conduction  band.  However,  this 
identification  was  questioned  by  Stratton  [52],  who  pointed  out  that  the 
peak  separation  for  bulk  states  emission  should  be  greater  than  the  gap 
energy,  and  the  energy  attribution  should  have  zero  amplitude  for 
energies  of  the  gap  region.  Russell  and  Litov  argued  that  the  poor 
resolution  of  their  analyzer  had  spread  the  expected  sharp  peaks  into 
the  broad  shapes  which  they  observed. 

Russell  and  Litov  [45]  apparently  encountered  some  experimental 
problems  [52].  These  authors  were  unable  to  measure  the  offset  of  their 
sample,  recording  only  the  transient  response  of  the  lock-in  amplifier 
on  an  oscilloscope.  Further,  they  claimed  to  see  no  symmetric  pattern, 
which  generally  indicates  either  contamination  or,  as  described  by 
D’Asaro  [14],  multiple  "tiplets."  No  information  on  vacuum  quality  is 
provided  in  their  paper.  The  FWHM  of  their  tungsten  distribution  was 
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reported  to  be  0.6 eV,  in  contrast  to  the  accepted  width  of  0.25eV  at 
300K. 

Rather  improved  studies  on  Si  were  made  by  Kleint  and  Kusch 
[36].  They  investigated  0.2  to  0.8  ohm- cm  n-type  Si,  of  (ill) 
orientation.  Using  a  high  quality  analyzer  designed  by  Young  [60], 
these  authors  described  three  distinct  cases.  First,  they  described 
seeing  a  single  peak,  showing  enlargement  of  width  and  energy 
displacement  with  increasing  field.  Secondly,  they  found  two  peaks, 
with  maxima  separated  by  about  1.2eV.  Thirdly,  they  reported  two  peaks 
separated  by  approximately  2.5eV.  Although  no  basis  was  given,  these 
authors  attributed  the  peak  of  lowest  energy  to  valence  band  emission, 
and  that  of  higher  energy  to  the  conduction  band,  as  had  Russell  and 
Litov. 


Previous  Work  on  Other  Semiconductors 

Work  on  GaAs  was  '■’erformed  by  Hughes  and  White  [30].  They 
used  heavily  doped,  n-type,  0.0001  ohm-cm  resistivity  material,  and  a 
field  desorbed  surface.  Using  a  Young  style  analyze.,  having  excellent 
resolution,  these  authors  found  a  large  wide  peak  located  almost  0.4-eV 
below  S  ,  and  a  small  subsidiary  peak  appearing  at  E  and  lower. 

Their  data  appeared  similar  to  this  author' s  own  data  on  heavily  n- 
doped  Si.  In  addition,  Hughes  and  White  described  small  sharp  peaks, 
located  between  E^  and  the  large  peak,  which  appeared  and  dis¬ 
appeared  on  consecutive  sweeps  of  their  retarder.  They  ascribed  the 
large  peak  to  the  valence  band,  and  the  subsidiary  peaks  to  surface 
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state  emission.  Those  authors  argued  that  conduction  band  emission  was 
absent  due  to  band  structure  effects,  even  though  GaAs  is  a  direct  gap 
material  and  transverse  momentum  considerations  would  be  of  negligible 
importance.  Further,  one  can  show  that  GaAs  of  the  quoted  resistivity 
implies  that  the  Fermi  level  is  within  the  conduction  band — the  sample 
is  degenerately  n-type. 

On  Ge,  several  researchers  have  published  results.  Arthur  [5] 

demonstrated  a  wide,  single  peak  appearing  below  E  .  This  peak  occa- 

Jr 

sionally  showed  small  oscillations,  especially  evident  at  low  tempera¬ 
ture.  Arthur  argued  that  such  a  wide  peak  was  consistent  with  a 
severely  degenerate  conduction  band  at  the  surface.  Also,  Arthur's 
results  seemed  independent  of  doping  type,  yielding  a  large  single  peak 
in  all  circumstances. 

Further  work  on  n-type  Ge  was  performed  by  Shepherd  and  Peria 
[47].  Attached  to  a  smaller,  wide  peak,  these  researchers  found  a  very 
strong,  sharp  peak,  of  width  0.1  eV.  The  sharp  peak  was  angularly 
localized  at  a  particular  bright  feature  of  the  annealed  Ge  pattern. 

For  other  emission  directions,  the  amplitude  of  the  subsidiary  peak 
varied.  The  strong  sharp  peak  was  ascribed  to  a  surface  state,  as  it 
does  not  fit  the  theory  for  the  Ge  (100)  direction  emission,  and  it  was 
sensitive  to  surface  contamination.  The  broader,  low  amplitude  peak  was 
argued  by  these  authors  to  be  valence  band  emission,  as  it  appeared 

below  E  and  had  the  expected  shape  of  valence  band  emission. 

2 

Lastly,  Kisker,  Mahan,  and  Reihl  studied  the  FEED  of  a  370i 
thick  semiconducting  EuS  film  deposited  upon  a  V  emitter  [34].  Their 
published  data  shows  a  large,  wide  peak,  located  nearly  3eV  below 
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E  ,  with  a  smaller,  subsidiary  peak  beginning  at  5  and  tailing 
Jr  F 

lower  in  energy.  These  authors  note  that  the  large  peak  shifts  toward 
lower  energy  and  broadens  with  increasing  field,  though  the 
subsidiary  peak  grows  in  relative  amplitude  and  remains  close  to 

bl. 

F 

Thus,  their  data  corresponds  in  general  shape  to  that  of  other 
researchers.  Kisker,  et  al.  argued  that  inelastic  scattering  of 
electrons  entering  the  EuS  layer  could  account  for  the  observed 
distribution.  Those  electrons  which  arrive  at  the  surface  unscattered 
by  traps  are  emitted  with  energy  equal  to  the  Fermi  level.  Inelastic 
scattering  reduces  the  energy  of  many  electrons,  however,  and  they  drop 
down  to  the  conduction  band  edge  or  to  trap  levels  within  the  conduction 
band.  Therefore,  a  large  electron  population  exists  at  and  near  the 
conduction  band  edge,  giving  the  large  peak.  As  the  applied  field  is 
increased,  they  argue,  the  large  peak  shifts  downward  with  the 
conduction  band,  while  the  small  peak  remains  near  the  Fermi  level. 
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Generalization  and  Discussion 


Comparing  all  of  the  accumulated  data  of  these  researchers, 

one  sees  a  general  trend.  Consistently,  a  single  large  peak  dominates 

the  energy  distributions.  This  single  peak  shifts  down  in  energy  with 

increasing  field,  and  a  subsidiary  peak  arises  just  below  E  .  With 

r 

increasing  field,  this  subsidiary  peak  increases  in  amplitude  compared 
to  the  lower  energy  peak.  The  subsidiary  invariably  appears  with  nearly 
pure,  high  resistivity  semiconductors,  or  with  heavily  doped  semicon¬ 
ductors  operated  at  low  (78K)  temperature.  This  implies  that  free 
carrier  density  is  strongly  related  to  the  PEED  characteristics. 

Although  the  explanation  put  forth  by  Kisker,  et  al.  [34]  is 
sensible  and  attractive  for  their  thin  film  data,  it  is  not  applicable 
to  the  other  samples.  For  example,  their  model  supposes  a  perfectly 
dielectric  semiconductor,  with  the  Permi  level  bending  down  parallel  to 
the  assumed  bending  of  the  valence  and  conduction  bands.  We  expect  band 
bending  of  the  order  of  pJx.  We  may  take  the  maximum  value  of  J  to  be 
1000  Amps/ cm2,  so  that  for  resistivities  p  of  order  .001  ohm-cm,  one 
volt  would  develop  for  an  emitter  of  length  one  centimenter.  This  is  an 
unreasonably  large  estimate  for  typical  emitters,  implying  that  this 
mechanism  is  unrealistic.  Por  this  reason,  we  must  discard  the  model  of 
Mahan,  et  al.  as  too  limited  to  explain  the  diverse  results  accumulated 
on  many  samples  of  high  conductivity  semiconductor. 
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Alternatively,  let  us  consider  another  model,  based  upon 
photoemission  work.  To  explain  their  photoemission  data,  Wagner  and 
Spicer  [55]  proposed  that  the  Si  surface  was  covered  by  a  large  density  of 
surface  acceptor  states.  For  their  cleaved  n+.  0.001  ohm-cm  'i  sample, 
they  concluded  that  the  density  of  surface  states  was  of  order  8  x 
10  /cm  —approximately  one  surface  state  per  surface  atom.  Moreover, 
they  concluded  that  Figure  39  is  an  acceptable  distribution  in  energy  of 
these  surface  acceptor  states.  From  that  figure,  we  find  the  surface 
states  overlapping  the  valence  band,  and  tapering  off  to  zero  density 
below  the  conduction  band  edge. 

This  picture  of  the  surface  state  density  leads  to  band 
bending  without  an  externally  applied  field.  Wagner  and  Spicer  cal¬ 
culated  that  the  band  bending,  for  the  planar  surface,  would  give  a 
0.2  eV  difference  from  the  flat- band  condition  at  a  distance  of  12  A 
from  the  surface.  This  bending  could  alternatively  arise  from  an 
imaginary  external  field  of  -0.51  V/A  [28].  Also,  there  is  no  net 
external  field,  due  to  the  surface  states,  because  the  band  bending 
causes  a  buildup  of  positive  charge  both  in  the  valence  band  and  on 
ionized  donor  states  near  the  conduction  band. 

The  application  of  an  external  field  will  flatten  the  bands. 
Assuming  that  the  surface  states  are  filled,  the  electrostatic  fields 
will  add  vectorially.  In  this  case  the  imaginary  internal  surface  field 
of  0.51  V/kA  will  add  to  the  emission  field  F/K,  with  K  the  dielectric 
constant.  The  band  curvature  will  be  V(x)  *  ’/^(FjexpC-x/L^).  Here, 

*  [  ekT/qaN0  ]  *  13  A.  The  field  just  within  the  surface  will  be 

dV(x)/dx  -  V0(F)/Ld.  Thus,  VQ(y)  •  (L  x  [  0.51V/A  -  ?]  /  K)  is  the 
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exteat  of  the  internal  band  bending.  As  shown  in  Figure  40(b),  the  exter¬ 
nal  field  of  order  0.3V/A  shifts  the  entire  surface  state  density  below 

E  ,  with  the  high  energy  tail  remaining  some  few  tenths  of  an  eV  below 

£ 

V 

This  analysis  of  Vagner  and  Spicer's  photoemission  data  sug¬ 
gests  an  alternative  mechanism  of  semiconductor  field  emission.  It 
suggests  that  emission  occurs  primarily  from  surface  states,  rather  than 
the  two  bands,  at  low  fields.  This  proposition  adequately  explains  why 
the  large,  major  peak  occurs  below  E  and  why  it  shifts  downward  in 

If 

energy  with  increasing  field.  It  explains  why  the  shift  is  not  attrib¬ 
utable  to  an  IR  drop  in  the  sample. 

Further,  the  invocation  of  surface  state  emission  may  explain 
the  discrepancy  in  shapes  of  the  FEEDs.  The  theoretically  predicted 
width  of  conduction  band  emission  is  very  narrow— narrower  than  the 
metallic  distributions  observed  for  V.  The  experimentally  observed 
distributions  are  always  much  wider  than  the  metallic  distributions. 

Also  the  theoretical  valence  band  distribution  does  not  change  width 
with  increasing  field;  it  only  changes  amplitude.  This  is  in  contrast 
to  the  experimental  observation  that  the  width  increases  with  increasing 
field.  In  addition,  the  wide  shape  of  the  experimental  distribution  may 
be  a  consequence  of  the  --creasing  density  of  surface  states  with  lower 
energy— similar  to  the  case  for  valence  band  emission,  which  also  gives 
a  wide  distribution.  Also,  a  change  in  the  surface  state  density, 
effected  by  annealing  or  desorbing  the  surface,  would  yield  a  change  in 
the  FEED.  This  is  observed  in  the  experimental  distributions. 

In  certain  circumstances,  emission  will  also  come  from  the 
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conduction  band.  This  will  occur  if  the  charge  on  the  surface  states  is 
insufficient  to  prevent  the  externally  applied  field  from  bending  the 
bands  downward  at  the  surface.  Such  will  be  the  case  if  emission  of 
electrons  from  the  surface  states  causes  their  population  to  diminish 
sufficiently.  Vagner  and  Spicer1 s  data  implies  that  the  bands  will 
never  bend  downward  for  typical  emission  fields  of  up  to  5x10  v/cm, 
in  the  absence  of  emission.  Indeed,  for  heavily  n-doped  semiconductors, 

we  do  not  generally  see  emission  at  energies  near  E .  However,  for 

£ 

lightly  doped  semiconductors,  or  heavily  doped  semiconductors  at  low 
temperatures,  we  do  see  emission  at  these  higher  energies.  Thus,  the 
concentration  of  electrons  in  the  surface  states  must  diminish  under 
these  emission  conditions,  to  about  0.8  of  the  original  surface  electron 
density. 

This  simple  model  employing  the  conduction  band  and  a  band  of 
surface  acceptor  states,  seems  to  explain  a  majority  of  the  experimental 
work.  The  experimental  reports  of  two  broad  peaks  can  be  ascribed  to 
these  two  sources.  This  model  easily  covers  the  observation  of  only  one 
peak,  that  of  the  lower  energy  surface  states,  as  well  as  that  of  the 
double  peaks.  It  also  handles  the  observed  fact  that  the  energy 
separation  between  peaks  increases  with  field,  and  that  the  higher 
energy  peak  increases  in  amplitude  with  increasing  field. 


CHAPTER  V 


PHOTO-FIELD  EMISSION 

Introduction 

The  effect  known  as  Photo-Field  Emission  (PFE)  or  alter¬ 
natively  as  Photo- assisted  Field  Emission,  was  found  in  Si  by  Busch  and 
Fischer  in  1963  [l 1 ].  It  had  been  discovered  much  earlier  by  Apker  and 
Taft  [4]  that  CdS  emitters  demonstrated  a  large  increase  in  FE  current, 
which  was  attributed  to  purely  to  photoconductivity  [4].  The  charac¬ 
teristic  non-linear  F-N  curves,  along  with  the  spectral  response  of  Si 
and  Ge,  implied  to  many  researchers  that  the  effect  was  also  totally 
attributable  to  photoconductivity.  Indeed,  typical  data  taken  from  a 
field  desorbed  Si  sample,  which  has  been  plotted  in  Figure  29,  shows 
that  the  effect  cuts  off  at  approximately  the  band  gap  energy. 

The  photoconductivity  argument  supposes  that  the  F-N  relation 
holds  true  at  the  surface.  The  F-N  equation  then  relates  the  total 

current  to  the  true  surface  potential,  V  .  V  is  related  to  the 

s  s 

applied  voltage  7^  by 

V  -  7  -  IR 

3  A 

R  -  R(T,A,I  ....) 

& 

Thus  the  F-N  equation  becomes 


(64a) 

(64b) 


I 

I 


I 
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I  ’  aF2  exp( -0.685  fP^/F)  (65) 

F  »  JS(V  -  IS)  -  / }V  (66) 

Solution  of  the  transcendental  equations  (65)  and  (66),  for  assumed 
typical  values  of  the  constants  a,  /J ,  R  yields  the  graphed  curves  shown 

1 

in  Figure  41.  Notice  that  the  shape  of  these  curves  closely  resembles 

the  experimentally  observed  graphs ,  shown  in  Figure  26,  when  the  results 

are  plotted  as  Log(l)  vs.  1/V  . 

A 

However,  the  results  do  not  appear  as  attractive  when  plotted 

in  F-N  fora  as  Log(l/V  )  vs.  1/V.,  as  in  Figure  42.  Although 

A  A 

the  data  appears  to  fit  theory  for  one  regime,  namely  low  V,  it  does  not 
fit  the  high  V  case.  Puther,  it  is  observed  that  p-type  Si  shows  a 
stronger  PFE  effect  than  n-type  Si  of  the  same  resistivity.  If 
photoconductivity  were  the  only  rationale  for  this  effect,  the  doping 
type  would  be  inconsequential.  In  addition,  researchers  found  that  the 
F-N  equation  is  not  obeyed  at  the  surface.  Using  a  retarding  potential 
analyzer,  the  surface  potential  can  be  found,  and  F-N  plots  of  the 
resulting  data  show  non-linearities.  Lastly,  the  PFE  current  always 
exceeds  the  dark  current,  even  with  identical  surface  potentials.  My 
data  on  this  latter  effect  is  shown  in  Table  2. 

t 

For  these  reasons,  the  photoconductivity  explanation  was 
rejected  for  a  more  rigorous  model  by  Yatsenko  [56].  Yatsenko  attempted 
to  model  FFE  as  an  extension  of  the  planar  MOS  structure  in  the  "quasi- 
dielectric"  limit.  This  planar  structure  admittedly  neglected 
geometrical  effects  of  the  actual  spherical  surface,  but  ris  felt  to 
involve  the  same  physical  properties  as  the  FFE  problem.  The  model 
presented  the  sample  in  terms  of  a  base  region,  a  depletion  region,  a 
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Figure  41 .  Log  of  Total  Current  as  a  Function  of  Inverse  of  Applied 
Voltage  for  Photoconductivity  Correction  to  F-N  Theory. 
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surface  inversion  layer,  and  a  surface  barrier  regime.  After  analyzing 
the  MOS  structure,  Yatsenko  concluded  that  under  all  typical  FE 
conditions,  the  internal  surface  field  would  be  higher  than  the 
avalanche  breakdown  field  in  Si,  typically  10^  V/cm.  Thi3  paradox  led 
Yatsenko  to  conclude  that  surface  states  must  therefore  be  screening  the 
sample  from  the  emission  field,  for  under  avalanche  breakdown  condi¬ 
tions,  he  reasoned  that  the  current  would  increase  unreasonably — 
behavior  not  experimentally  observed. 

The  conjecture  of  surface  states  screening  the  sample  is 
reasonable.  As  demonstrated  by  my  data,  presented  in  the  section  on 
FEED  of  Si,  the  energy  lowering  of  the  main  peak  with  increasing  field 
conforms  with  this  interpretation.  Further,  the  width  of  the  distri¬ 
bution,  even  for  heavily  doped  Si,  may  indicate  that  an  increasing 
density  of  surface  states  with  lower  energy,  as  discussed  in  that 
section.  In  addition,  Stetsenko,  Yatsenko,  and  Miroshnichenko  [50] 
found  that  a  single  applied  voltage  pulse  broke  into  two  current 
segments  of  discernably  different  time  constants;  presumably  the  two 
were  due  to  emission  from  bulk  states  and  from  surface  states.  Lastly, 
the  two  peak  amplitudes  showed  hysteresis,  interpreted  to  be  further 
evidence  of  the  surface  states  having  a  larger  time  constant  than  bulk 
states. 

In  spite  of  the  reasonable  argument  that  surface  states  are 
screening  the  bulk  from  the  external  field,  Yatsenko  overlooks  several 
important  points.  First,  the  observation  of  PFE  coincides  with  a  shift 
in  the  dark  surface  potential.  Second,  the  large  sample  resistance 
limits  the  current,  even  if  avalanche  breakdown  occurs  in  the  surface 
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region.  Both  these  points  suggest  that  the  sample  resistance  must  be 
taken  into  account  more  explicitly  than  in  Yatsenko's  treatment. 

Further,  FEED  data  on  low  temperature  Si,  and  the  hysteresis  data  of 
Stetsenko,  et  al.,  indicate  that  the  surface  states  deplete  to  some 
extent  during  PFE  emission. 

This  last  point  indicates  that  Yatsenko  was  correct  in  his 
initial  assumptions  about  the  mechanisms  of  PFE,  namely  that  PFE  occurs 
at  the  point  where  the  surface  states  become  depleted.  However,  his 
theory  showed  that  avalanche  breakdown  must  occur  at  the  surface  region 
if  the  surface  states  deplete  substantially.  Yatsenko  reasoned  that 
such  a  breakdown  must  necessarily  be  catastrophic,  as  current  is  a 
continuous  quantity  throughout  the  sample.  Plagued  by  the  conclusion  of 
avalanche  breakdown,  he  was  unable  to  forfeit  the  role  of  surface  states 
in  the  model.  In  other  words,  if  the  surface  states  alone  control  the 
emission  characteristics,  then  one  must  be  left  with  the  pure  photo¬ 
conductivity  argument.  However,  if  both  surface  states  and  the  bands 
contribute  to  emission,  then  additional  current  terms  come  into  play. 


A  Second  Model 


Rather  than  revise  Yatsenko's  model,  I  propose  the  following, 
second  model.  This  model  simplifies  the  PFE  by  considering  emission 
from  two  energy  regimes.  In  reality,  three  energy  regimes  are  possible — 
the  surface  states,  and  two  bulk  bands.  However,  as  no  knowledge  is 
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presumed  of  the  actual  form  of  the  surface  state  wavefunctions,  I  will 
take  the  valence  band  as  a  "reasonable"  substitute  for  the  surface  state 
emission.  Further,  I  will  assume  that  the  conduction  electrons  may 
possess  a  zero  minimum  transverse  energy,  i.e.  a  free  electron, 
spherical  energy  surface  for  the  conduction  band. 


Account  of  Photoconductivity 

To  take  account  of  photoconductivity,  we  will  again  assume 
that  the  surface  potential  V  and  the  field  F  are  directly  proportional. 

3 

Thus, 

F  -  yS(V  -  IR.)  -  /}V  (67) 

A  3 

and 

s  1  /  (5  r  )  (68) 

v 

In  determining  R  ,  we  use 

-  R (  sample  dimensions)  (69) 

P\  -  a\  M  °o  *  (70 ) 

a\  *  (nem„  +  pe^p  )  +  (dneftn  ♦  Apejup  )  (71  ) 

now  An  -  Jp  -  A1 (Intensity  I  ,  A .Reflection  coeff. , 

n 

a  (A)...)  (72) 

n-p  -  nt2  -  1.4  x  l020/cm6  x  (T/300K)3  (73) 

For  p-type  Si,  p  ■  H  ,  n  -  n2 /R A«  At  78K,  n  -  3.44  x 
1018/R 

A 

For  lightly  doped  Si,  say  1014/cm5,  this  gives  n  -  3.44  x 
104/cm^.  So 
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A _ gdark  =  ne<i"  +  Pe^P 

^dark  °X  (n+A’)e«n+  (p+A‘)e^p  (74a) 

*  H^eup/fA  ' e,an  +  (N^  +  A')e,up]  (74b) 

with 

"dark’  '/"dark  (75) 

To  determine  A  * (i  ,.l,R,...),  we  assume  that  each  generated 
n 

carrier  has  a  lifetime  t.  The  reflection  coefficient,  R,  is  taken  to  be 

0.3.  Complete  absorption  is  assumed— in  reality  absorption  is  related 

to  the  sample  thickness  and  the  absorption  coefficient  at  the  wavelength 

of  the  incident  radiation.  With  these  approximations, 

A*  *  (#e“/cin)  *  (#e  generated/ cm  -sec)  x 

(2/-r  rL  /tt2L  )x(1-R)xt  (76) 

o  0 

A'  -  [i  A  5.03  x  1014/cm2-sec]  x  [2/r]  x0.7r  (77) 

with  A  in  A,  Iq  -  intensity  in  Watts/cm  .  The  variable  "r  is  an  average 

radius  of  the  sample  *  (r  +  r.)/2,  in  cm.  In  this  approximation, 

b  t 

effects  of  surface  absorption,  bulk  absorption,  etc.  are  ignored.  We 
aim  here  only  to  arrive  at  a  functional  approximation  to  the  absorption. 


Note  on  the  Lifetime 

The  recombination  of  csrriers  follows  the  law 

dn/dt  ■  -Cn«  p  (78) 

with  "  a  proportionality  constant.  If  p  *  N  ,  the  solution  is 

a  -  n  exp(-CN  t)  -  n  exp(-t/r  ) 

O  A  0 


(79) 
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30  T  “  Tn  a  l/CN  .  This  means  that  the  reduction  in  carrier  lifetime 
A 

eliminates  or  substantially  reduces  the  PFE  effect  for  heavily  p- 
doped  samples. 


Sample  Resistance 


We  determine  B  M  from  the  classical  expression  8  *  pL/A 
applied  to  a  cone  shaped  sample  of  length  and  having  a  base  and  tip 
of  radii  r^  and  rt  respectively.  The  angle  of  the  aone  is  Ym  (rfc  - 
rt)/LQ.  Applying  the  resistance  formula  yields 

R.l  ’  fyJo*dL/[A(L)]  -  (p\/n)  ./^dL/t  (rb  -  VL)2]  (80a) 

*  (px/nY)  x  (1/r  -  1/r  ) 

X  D 

“  P\ L  /nr  r. 

•x  o  t  b 

or  *  (r^>>rt)  (80b) 

The  contribution  of  the  hemispherical  cap  on  the  resistance  is 
negligible,  as  the  path  length  involved  is  only  of  order  r  <<  L  . 


Photogeneration  of  Carriers 

The  photogeneration  of  free  carriers  does  more,  qualitatively, 
than  just  change  the  sample  resistivity.  Free  electrons  are  accelerated 
toward  the  surface  by  the  internal  field.  These  photogenerated 
electrons  have  higher  total  energy  than  those  in  surface  states  or  the 
valence  band.  Thus,  these  conduction  electrons  may  be  emitted 


102 


preferentially,  possessing  substantially  higher  tunnelling  probability. 
If  not  emitted,  these  electrons  may  become  trapped  at  the  surface, 
cascading  into  the  lower  energy  surface  states. 

We  proceed  along  these  lines  as  follows. 

Let 


total 


Jdark(A)  *  Jlight(A) 


(SI  ) 


where  J ,  , (.1)  is  the  valence  band  current,  calculated  with  Ri  used, 

dark 

Jlight(A)  is  the  contribution  from  the  electrons  generated  by  light.  We 
approximate  the  dark  current  by  the  valence  band  contribution, 
neglecting  surface  state  emission,  taking 


J  (A)  -  J  (F) 

dark  v 

(82) 

and 

I  *  J  ^  ,  x  Area  of  tip. 
total 

(83) 

In  turn,  J^f)  is  evaluated  using  the  Stratton  model  and  free  electron 

mass  for  the  holes.  This  calculation  may  be  approximated  by 

Jy(F)  -  8*mef(E)  3«eF2/[9h5(2m)1 /2  4  (X  +  SQ)]  x 

exp[  -4  (2  m) 1  /2  ( X  +  E,  9  /2  /5UeF]  (84a) 

or,  Jy(F)  “  [7.6  x  10^Amps/cm2]  x  F2  x 

(l/(.V  +  Ea))  x  exp[-0.685(X+  E9)3/2/f]  (84b) 

with  X,  Eg ,  in  eV,  and  F  in  V/A. 

To  evaluate  J  (A),  we  assume  that  all  electrons  are  at  the 
light  ' 

bottom  of  the  conduction  band.  Therefore,  we  may  express 

Light  •  *"• °f  1  'W-0 

-Cl -R]  x  [# photons  generated/ cm2 -sec]  x 
[Active  area]  x 
[Probability  of  emission]  x 

[l.6  x  10_1  ^ coui/ photon  absorbed]  (85) 
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The  components  of  (85)  are  evaluated  below. 

(l  )  R  *  reflection  coefficient.  Assuming  normal  incidence  of  un- 

2 

polarized  light,  R  =  [ ( ( n2  -  u1  )  +  k2)/((n2  +  n1  )2  +  k2)]  .  The  ex¬ 

tinction  coefficient,  k,  is  negligibly  small  for  photons  of  energy  below 
1.2eV.  Here,  n1  -  1  for  vacuum,  n2  “\/*”  3  3.435  for  Si.  So  ( 1  — R )  » 
4n1n2/(n1  +  n2  )  -  0.7. 

(2)  #Photons  generated/cm2- sec.  This  is  the  intensity  I  divided  by  the 

energy  per  photon,  or 

[#photons/cm2-sec]  =*  I  A  (A)/hc 

n 

*[5.03  x  10UA(1)  Iq( Watts/ cm2)]/(cm2-sec) 

(3)  Active  area.  This  is  the  area  of  the  sample  in  which  generated 
electrons  reach  the  surface.  If  we  assume  a  recombination  lifetime,  r, 
for  the  excess  electrons,  then 

n(t)  -  n0exp(-tA) 

is  the  number  left  at  time  t  after  generation  at  time  t»0.  Assuming  a 
drift  velocity  vd  -  /iE,  we  get  x  ■  vdt  and 
n(x)  -  n0exp(-x/vdT)  . 

The  number  arriving  at  the  surface  x*0  from  a  distance  x  away  will  be  N# 

CO 

*  / n(x)W(x)dx.  Here,  W(x)  ia  the  width  of  the  conical  shaped  tip  of 
angle  defined  just  before  Equation  80.  Thus,  W(x)  *  (2rt  *  2  Vx). 
Integration  yields 

ff  -  2n0vdT(r,  +  Fvdr). 

Therefore,  the  effective  area  is  JIg/n8,  or  [Active  area]  -  2vdT(rt  x 
vdKf)  or  (rb  x  L0)/2,  whichever  is  less. 

(4)  Probability  of  emission.  We  will  take  this  to  be  the  WKB  solution 
at  the  conduction  band  energy,  neglecting  the  image  potential  and  band 
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structure  effects,  and  assuming  the  free  electron  solution.  Thus, 
[Probability  of  emission]  •  D(EC  )  *  exp(-0. 6a5J3/^F) . 

(5)  [Surface  emission  area].  We  take  this  to  be  2ttt  2 ,  neglecting  ob- 
vious  structure  characteristics,  such  as  facetting,  which  reduces  the 
emitting  area  of  the  surface.  The  field  desorbed  surface  appears  to  have 
uniform  emission,  and  is  more  closely  represented  by  this  expression. 

Therefore,  substituting  the  above  into  equation  (85)  yields: 

Jlight(A)  "  ^aA(1,35  x  10'4)vdT]  x 

r  O 

[r  +  VvdT]Amps/cin  x 

(exp[-0.685  X3/2/F])/[2«  r  2]  (86) 

or  *  [l  A  (1.35  x  10“4)/4]  x  [r.  x  L  ]  x 

n  do 

exp[ -0.685 -V3/%]/[2  n  r  2] 

X 

To  this  point,  we  have  a  formulation  of  PFE  embodied  in  equations  (67) 
through  (86).  Graphs  of  data  appear  in  Figure  43,  with  the  approxi¬ 
mations  taken  below: 

I  ■  I  x  1  Watt/ cm2  E  ■  1.1 eV 

n  n  g 

Ha  -  10U/cm3  A  -  104A 

“  lO^cm/sec  X  *  4eV 

T  “  10“^aec  x  (10^/tt,)  L  *  5mm  “  0.5cm 

A  o 

r  -  3000A  -  3  x  10~5cm  ,“p  -  500cm2 /V-sec 

X 

r  *  5  x  10“3cm  P n  *  1500cm2 /V-sec 

D 

V  “  (r  -  r.)/L  *  10  radians 

b  to 
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With  these  values,  the  set  of  equations  reduce  to  the  five 

below: 

y(v/A)  *  (V  I  103  -  I  I  RA)/(1.5  X  io4i) 

A 

Jlight^A/,Cm2 ^  “  0*245  x  10^)  x  1^  x  exp[-5.48/?] 

Jdark(A/cn,2)  “  (l*49  x  1C>9  x  ?2)  x  «pC-7-89/p] 

lU>  '  '•’light  *  '’dark’  1  <5.65  ,  10‘9) 

RA(ohms)  -  (1.33  x  108)  x  [l  +  t.13  x  I  ]_1  (87) 

n 

Ve  observe  from  Figure  43  that,  while  this  formulation  of  the  model  does 
give  a  PFE  current  with  light  incident  that  exceeds  the  dark  current  at 
the  same  surface  potential,  it  does  not  yield  the  up- turn  behavior  at 
high  V  near  the  origin  of  the  abscissa.  For  this  reason,  we  must 
postulate  that  a  breakdown  mechanism  exists  at  high  cross-sample 
voltages. 

Two  types  of  breakdown  mechanism  seem  possible.  The  first, 
avalanche  breakdown,  may  occur  when  the  internal  field  within  the  Si 
reaches  the  breakdown  field.  For  Si,  this  breakdown  field  is  a  function 
of  doping,  and  is  of  order  lO^Y/cm  for  a  background  doping  of 
1014  /cnP .  Avalanche  breakdown  will  multiply  the  expected  current  at  a 
particular  cross-sample  voltage  by  some  factor  M.  This  implies  that  the 
cross-sample  resistance  Rj  should  be  multiplied  by  the  factor  1  /M. 
However,  the  avalanche  breakdown  mechanism  reduces  the  resistance  only 
in  the  region  of  the  sample  which  is  actually  supporting  the  breakdown 
field — the  rest  of  the  sample  must  retain  the  same  resistivity. 
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Therefore,  we  may  envision  that  the  surface  region  of  the  sample  will 
breakdown  earliest,  and  that  a  breakdown  region  will  grow  in  length  back 
toward  the  sample  base. 

The  second  mechanism  is  punch- through.  Punch- through  is 
typically  a  bipolar  transistor  breakdown  mechanism,  in  which  the  base 
region  is  covered  by  the  depletion  regions  of  the  emitter  and  collector. 

In  our  crae,  the  depletion  region  extends  from  the  surface  back  toward 
the  base  of  the  sample,  as  described  above. 

Slight  differences  exist  in  the  application  of  these  two 
mechanisms.  In  the  second,  the  depletion  region  is  assumed  to  have  zero 
resistance.  In  avalanche,  however,  the  resistivity  is  merely  reduced  to 
1/M  of  its  original  value.  The  avalanche  mechanism  will  be  utilized  in 
the  breakdown  mechanism  in  this  model. 

Dielectric  Model 

The  simplest  model  -of  the  sample  is  that  of  a  pure  dielectric. 

Thus  we  will  get  a  contribution  to  the  dark  current  from  the  avalanche 

J 

throughout  the  sample. 

To  compute  the  avalanche  current,  we  assume  that  the  non- 
avalanched  current  is  multiplied  by  some  factor  M.  As  this  occurs  at 
some  fixed  cross-sample  voltage  V,  the  original  sample  resistance  R  must 
be  reduced  to  a  value  R/M.  Further,  some  of  the  avalanche  generated 
electrons  go  into  the  conduction  band  where  they  may  be  emitted. 

To  evaluate  the  number  entering  the  conduction  band,  we  note 
that  the  conductivity  increases  by  the  same  factor  M,  so 

erM  *  <7  +  (M-1  )o  ■  a  *  (excess  c r) 

i 
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*  a  *  ne/in  ♦  Pe^p 

with  n  and  p  representing  excess  carriers  generated  by  the  avalanche 
process.  We  assume  that  *  3jUp,  and  that  the  number  of  excess  holes 
generated  equals  the  number  of  electrons  generated.  Thus, 

(M-1  )o  «  4pea 


or  n  *  p  =*  (M-1  W4e|Up 

Using  <r  -  [if  ♦  A']e/iP+  [a  '  ♦  n±  /N  ]eu„ 

Combining  terms, 

a  *  [»A  ♦  +  4A*]e«p 

Thus  n  -  (M-1  )/4  x  [S  +  3n.2/N.  +  4A'] 

A  X  A 

The  number  of  these  n  electrons  which  arrive  at  the  surface 


per  unit  time  per  unit  area,  N,  must  be  evaluated.  Using  a  Maxwellian 


distribution  of  velocities  leads  tos 


With  A  -  m/2kT, 

^(v)  *  ^A/*]^2  x  exp(-A-v2] 

to 

Thus,  H  -  n[A/«]  x  Jv^exv(-Av^)  dvx  x 

0 

/ex p(-Avy2)  dv  x  /expC-Av^2  )  dv^ 

-•  -O 

so  N  -  n[A/*]5/2[l/2A]|>/A] 

yielding  H  *  nfkT^^m]1^2 

We  also  require  an  expression  for  M  itself.  Bar-Lev  gives  an 
empirical  expression  for  diode  breakdown  as  [7]: 

M  -  1  /  [l  -  (S/E^)2®] 
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The  variable  m  is  empirically  determined,  and  is  usually  between  3  and 
6.  To  a  first  approximation,  we  may  take  the  internal  field  to  be  the 
total  IR  drop  across  the  sample  divided  by  the  total  length  of  about 
5mm.  Then,  taking  for  Si  of  light  doping  to  be  2  x  105V/cm,  and 
nr*3,  we  get: 

M  -  1  /  [l  -  (IR/ (0.5cm  x  2  x  105V/cm))6] 
or  M  »  1  /  [l  -  (IR/105)6] 

The  effect  of  choosing  a  multiplier  function  which  is  less  "hard,"  such 
as  M  *  [l  -  IR/10^]“^  will  also  be  demonstrated.  The  model  employing 
the  latter  equation  will  be  termed  "soft"  breakdown. 

Thus,  the  additional  current  density  expected  from  the 
avalanche- generated  conduction  band  electrons  will  be 

J  «  n  x  [kT^Trm]1  ^  x  [l.6  x  10  ^coul/chg.]  x 

IwD 

[probability  of  emission]. 


Pull  Model  with  Breakdown 

The  full  model,  with  breakdown  now  included,  is  embodied  in 
the  following  set  of  equations.  Here,  the  previous  parameter  assump¬ 
tions  have  been  inserted. 

(l  }  Jdark  “  1  x  109  x  F2  x  exp(-7.89/F) 

(2)  Z  -  ([1  .245  x  105  x  I  x  M']  ♦  [a.  32  x  10~13  n] ) 

light  n 

x  exp(-5.48/F) 
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(3) 

n  =  [  (M-1  )/4]  x  [10U  ♦  l 

(4) 

A '  -  [2.82  x  1013  x  I  1 

n 

(5) 

F  =  [V  x  103  -  T  x 

A  total 

(6) 

Ra=  [r0/m]  /  [l  +  1.13  X 

(7) 

(Jj  ,  +  J,.  J  x  5.65  : 
dark  light 

(8a) 

M=  M.  =  [l  -  (2  x  I 

1  total 

(8b) 

M»  M,  -  [l  -  (2  x  I 

2  total 

(9) 

M*  »  [n  ♦  A*]  /  A’ 

(10) 

r  =  3 000 A 

/  5rt 


'total 


RA  /  (2  x  105) )6 ]_1  "hard” 


x  Ra  /  (2  x  105))  ]'6  "soft" 


The  model  is  most  easily  solved  in  the  following  order.  First, 
one  chooses  values  for  the  intensity,  I  ,  dark  resistance  R„,  and 
arbitrary  total  current  value  Secondly,  equations  (6)  together 

with  either  (8a)  or  (8b)  are  solved  iteratively,  yielding  RA.  Using  RA, 
M  is  obtained  from  (8a)  or  (8b)  appropriately.  Lastly,  the  remaining 
equations  are  solved  iteratively,  using  equation  (7)  as  the  test. 

Results  of  this  model  appear  in  Figures  44  and  45-  R0  has 
been  taken  to  be  lO1^  ohms  in  both  cases.  Figure  44  shows  the  result  of 
using  ,  the  "hard"  breakdown  expression.  In  contrast,  Figure  45  shows 
the  "soft"  breakdown  case,  utilizing  M^.  The  hard  breakdown  quali- 
tatively  follows  the  experimental  data,  showing  ln(l/V  )  with  a 

A 

lcoal  minimum.  Alternatively,  the  "soft"  breakdown  gives  the  smooth 
transition  between  high  and  low  current  regimes,  apparent  in  the 
experimental  data,  compared  to  the  sharp  slope  of  the  "hard"  case  at  low 

(1/V.)  values.  Some  combination  of  the  two  approaches  used,  or  a  better 

A 

empirical  expression  for  X  would  no  doubt  yield  better  qualitative 
agreement  to  the  empirical  data. 
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Figure  44.  Fowler-Mordheim  Plot  of  Revised  Theory  with  Hard 
Breakdown  Mechanism. 


•  -***»*» 


113 


Lastly,  it  should  he  noted  that  the  graphs  derived  from  these 
models  show  discrepancy  in  the  low  current  regime.  Experimentally,  the 
FEED  studies  show  emission  generally  originates  from  one  wide  energy 
regime,  ascribed  to  surface  states.  For  high  current  densities,  a 
subsidiary  peak  at  higher  energy  appears,  but  it  is  of  lower  intensity 
than  that  of  the  surface  state  emission.  In  these  models,  the  author 
has  neglected  treatment  of  the  surface  barrier,  present  with  negatively 
charged  surface  acceptor  states.  For  low  current  densities,  it  is 
envisioned  that  this  barrier  will  substantially  reduce  the  conduction 
band  contribution  to  the  current,  and  eliminate  the  discrepancy  present 
in  these  models. 

Taking  this  into  account,  theoretical  plots  using  the  same 

model  but  with  J,  ^  reduced  by  a  factor  of  5GOO,  are  shown  in 
light 

Figures  46  and  47.  These  figures  indicate  that  the  discrepancy  at  low 
current  levels  is  eliminated  by  neglecting  the  conduction  band 


contribution 


in  KV) 
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Figure  46.  Fowler-No rdheim  Plot  of  Revised  Theory  with  Hard 

Breakdown  Mechanism  Without  Conduction  Band  Contribution. 
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CHAPTER  VI 


CONCLUSION 


la  summary,  the  thesis  has  examined  the  results  of  both  past 
and  current  experimental  results  with  respect  to  the  present  theory  of 
FE  from  semiconductors.  In  general,  application  of  the  Stratton  theory 
predicts  either  one  or  two  peaks  in  the  FEED,  corresponding  to  the 
valence  and  conduction  bands.  The  separation  of  these  two  peaks  at  the 
base  is  taken  to  be  a  constant,  the  thermal  gap  energy,  if  surface 
states  are  neglected.  Further,  as  the  Fermi  level  is  shifted,  the 
valence  band  peak  theoretically  should  not  change  shape — only  amplitude. 
Experimentally,  both  the  peak  separation  and  lower  energy  peak  shape 
widen  with  increasing  field.  The  single  emission  peak,  seen  for 
degenerately  n-doped  semiconductors  at  300K,  always  appears  some  tenths 
of  an  eV  below  the  Fermi  level,  even  at  very  low  current  densities.  A 
smaller,  subsidiary  peak  appears  at  78K  for  these  heavily  n-doped 
samples  at  high  fields.  The  subsidiary  peak  also  appears  in  data  of 
other  researchers  for  lightly  doped  semiconductors,  and  those  at  low 
temperature. 

Comparison  with  photoemission  results  leads  to  a  different 
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model  of  PE  from  semiconductors.  The  data  of  Wagner  and  Spicer  [55]  on 
photoemission  from  heavily  n-doped  Si  shows  a  wide  band  of  surface 
states  contiguous  to  the  valence  band  and  tailing  off  to  zero  a  few 
tenths  of  an  eV  below  the  conduction  band  edge.  Combining  thi3  inter¬ 
pretation  with  the  application  of  an  electron  emission  field  leads  to 
the  conclusion  that  surface  state  emission,  rather  than  valence  band 
emission,  is  responsible  for  the  large  peak  observed  below  the  Fermi 
level  on  heavily  n-doped  semiconductors.  This  conclusion,  implying 
surface  state  emission,  explains  the  energy  position  of  the  lower  peak 
and,  qualitatively,  its  width.  The  upper,  subsidiary,  peak  is  tenta¬ 
tively  ascribed  to  emission  of  conduction  band  electrons. 

Such  a  conclusion  from  FEED  studies  and  photoemission  also 
concurs  with  work  on  PFE.  For  example,  Stetaenko,  fatsenko,  and 
Miroschnichenko  [50]  observed  two  separate  time  constants  from  single 
voltage  pulses  applied  to  emitters.  In  addition,  they  found  hysteresis 
in  the  emission  characteristics  of  pulse  pairs;  the  second  of  the  two 
pulses  always  showed  reduced  emission.  To  these  authors,  their  results 
indicated  that  emission  occurs  primarily  from  surface  states,  which 
deplete  under  strong  emission. 

Surface  states  may  also  play  a  role  in  explaining  anomalies  in 
FE  work  function  determinations.  Although  this  point  was  not  inves¬ 
tigated  in  the  thesis,  Busch  and  Fischer  [ll]  determined  that  work 
function  values  seemed  to  correlate  to  the  Fermi  level  in  the  bulk. 

That  is,  an  n-type  emitter,  having  a  high  bulk  Fermi  level,  also 
demonstrated  a  large  work  function.  This  thesis  presents  data 
indicating  that  field  desorbed  and  annealed  surfaces  show  different  work 
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function  values.  Perry  [42,43]  found  that  Si  emitters  gave  two  or  more 

distinct  work  function  values  for  the  same  emitter.  Generally,  the  » 

change  in  work  function  came  at  high  current  densities.  His  data  shows 

that  the  value  of  work  function  may  increase  or  decrease  upon  changing 

into  the  high  current  regime.  Also,  Allen  [ 1 , 2 ]  observed  that  a 

hysteresis  existed  for  some  of  his  t normally  annealed  Si  emitters.  The 

surface  state  density  may  thus  be  changing  in  these  various 

circumstances. 

In  application  to  a  theory  of  PFE,  Yatsenko  [56]  was  the  first 
to  take  account  of  surface  states.  His  theory  initially  assumed  that 
such  states  depleted,  allowing  field  penetration  into  the  semiconductor 
surface.  However,  he  abandoned  this  notion  after  showing  that  the 
penetrating  field  was  greater  than  that  necessary  for  avalanche 
breakdown.  In  Yatsenko' s  model,  avalanche  necessarily  leads  to  an 
infinite  current.  The  experimental  FEED  work  indicates  that  some  field 
penetration  does  exist  for  low  temperatures  and  dopings— precisely  the 
case  for  FEE.  Thus,  a  new  model,  eliminating  or  reducing  the  breakdown 
problem,  was  needed  to  link  the  FEED  evidence  to  the  PFE  results. 

Such  a  model  is  suggested  in  this  thesis.  In  this  model,  the 
photon  flux  generates  carriers  with  two  functions.  The  carriers  reduce 
the  sample  resistance,  and  electrons  generated  increase  the  conduction 
band  population.  Although  breakdown  is  seen  to  be  necessary  to  give  the 
correct  shape  to  the  PFE  F-N  plots,  the  breakdown  does  not  yield 
infinite  current.  This  latter  is  because  the  sample  resistance  never 
drops  to  zero.  Thus,  the  surface  field  is  limited  by  the  potential  drop 
within  the  sample.  As  a  result,  the  theoretical  F-N  plot  shows 


qualitatively  the  same  behavior  as  that  of  the  experiment.  It  is 
suggested  that  a  better  correspondence  between  theory  and  experiment 
would  result  from  a  more  sophisticated  treatment  of  the  breakdown 
mechanism. 
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APPENDIX  A 


3ASIC  PROGRAM  FOR  EXPERIMENTATION 

The  Basic  program,  contained  in  this  appendix,  controls  the 
FEED  experiment.  After  loading  this  program  into  the  Apple  II,  the  user 
must  next  load  in  the  machine  language  program  contained  in  Appen-dix  B. 
Typing  ]RUN  causes  a  "menu"  of  options  to  appear  for  the  program 
user. 

These  options  appear  in  lines  130  to  202.  After  a  choice  is 
made  by  the  user  in  line  210,  control  is  transferred  to  specific  program 
sections,  as  defined  in  lines  215  to  274.  Such  sections  prompt  the  user 
for  input  whenever  decisions  are  required.  Musical  prompts  are  incor¬ 
porated  to  direct  user  attention  to  equipment  problems  or  when  physical 
manipulation  is  required,  such  as  during  the  X-Y  plotting  routine. 

Once  data  is  taken,  it  can  be  stored  on  cassette  tape.  The 
user  first  must  determine  the  extent  of  memory  which  the  data  occupies. 
This  is  accomplished  by  typing  "D"  from  the  menu.  The  resulting  prompt 
character,  a  star  (*),  shows  that  the  computer  is  in  the  "monitor”  mode. 
Typing  *7000  next,  then  holding  the  REPEAT  key  down  with  the  RETURN  key 
depressed  simultaneously  causes  the  memory  locations  to  appear  rapidly 
on  the  screen.  If  they  are  zeros  until  *8000  is  reached,  then  the 
minimum  value  of  the  data  should  be  *8000.  Otherwise,  the  minimum 
should  be  *7000.  The  user  determines  the  maximum  value  by  testing 
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memory  values  higher  than  *8000.  Typically,  one  teat  *9000,  *A000,  etc. 
until  zeros  are  encountered.  The  end  of  the  data  is  determined  by 
testing  between  the  known  zero  and  known  non-zero  locations.  The  user 
starts  the  cassette  recorder  in  the  "Record"  mode,  then  types 
•mini.maxiW  to  Write  onto  the  tape  from  minimum  to  maximum  locations. 

At  this  point,  the  user  should  make  notation  of  the  stored  memory 
minimum  and  maximum.  To  reload  data  stored  on  tape,  the  user  types  "D" 
from  the  menu.  After  starting  the  cassette  recorder  on  "Play,”  the  user 
types  *8000.97P0R,  for  example,  to  Ilead  the  recorded  data  into  locations 
*8000  to  *9700.  It  is  important  to  realize  that  the  data  must  be  placed 
back  into  the  same  memory  locations  from  which  it  was  originally  stored. 
Also,  it  should  be  noted  that  the  cassette  recorder  should  have  its 
volume  and  tone  controls  both  set  to  the  value  9- 
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jlisi 


5  GOTO  180 

6  FOR  I  *  23672  TO  3  27  €.7 :  POKE  I ..  0s  NEXT  Is  GOTO  1 00 
10  INPUT  "NEW  DATE  OP  SAMPLE?"! A*:  IF  A*  <  >  "V"  THEN 

35 

15  INPUT  "DATE* • M0,DV,7P>"!M0,DV,VR:  POKE  27484,  MO*  POKE 
27485 , DV i  POKE  27486 . VP 

21  PRINT  "SAMPLE  TYPES* 1  =  W":  PRINT  TAB<  13>"2  »  SI 
LI CON.  INTRINSIC"!  PRINT  T AB <  13>"3  »  SI  P-TVPE” 
s  PRINT  TAB'-'  1  ’  '  "4  a  SI  N— TYPE" :  PRINT  TAB<  13> 

"5  =  OTHER  METAL  :  PRINT  T AB 13>"6  =  OTHER  SEMI 
CONDUCTOR" 

25  INPUT  "WHICH  SAMPLE  TYPE?": 31*  POKE  27483,51:  INPUT 
"ON  GRAPHING-WHhT  MESH  SIZE?"! A? 

35  POKE  32779, <.  PEEK  <27404>>:  POKE  32738, <  PEEK  <274 

05  >  >*  POKE  32731, <  PEEK  < 27486 >  }  :  POKE  32773, <  PEEK 
< 27483 >.>*  POKE  6,  AS!  RETURN 
59  A1  a  0 

68  POKE  -  16304,8!  POKE  -  16299,0:  POKE  -  16297,0* 
POKE  -  16301,8 

70  INPUT  "WANT  OUT?"!Al»*  IF  A2  *  0  THEN  A1  =  8 
72  I F  A2  —  1  THEN  A1  a  1 
75  IF  Al*  a  "V"  THEN  GOTO  100 

38  IF  Hi  =  1  THEN  POKE  -  16380,0:  POKE  -  16381,8 

35  IF  A1  =  0  THEN  POKE  -  16299,0:  POKE  -  16382,0 

37  IF  Hi  =  0  THEN  A2  «  1 

33  IF  Al  =  1  THEN  A2  *  0 
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GOTO  70 

100 

TEXT 

128 

PRINT  " 

MENU: TYPE  ’F’  FOR  FIXED  RETARDER  UALUE 

":  PRINT 

:  PRINT 

TAB<  11J"’ I'  FOR  INFO  ON  PARAMETERS" 

:  PRINT 

130 

PRINT 

TAB<  1 1  >  "  ’  M  ’  FOR  PUTTING  RETARDER  AT": 

PRINT 

TAB<  17> "  10  'JOLTS":  PRINT  TAB<  11>"’G’  FOR 

RUNN 

ING  THE 

EXPERIMENT" 

148 

PRINT 

TAB<  ilV’U’  FOR  RAW  NUMERICAL  DATA": 

PRINT 

198 

PRINT 

TAB<  1 1 > " ' L '  FOR  LOOKING  AT  A  GRAPH": 

PRINT 

TAB<  1 

7> "OF  THE  DATA" 

282 

PRINT 

TAB<  1 1 > ” ’ D ”  FOR  READING-WRITING  DATA" 

:  PRINT 

:  PRINT 

TAB<  ll.V’C’  FOR  COMPARING  DATA":  PRINT 

218 

PRINT  " 

" :  INPUT  "WHAT  COURSE?" ! Zl* 

215 

IF  Zl* 

=  "I"  THEN  1350 

22© 

IF  Zl* 

=  "F"  THEN  290 

230 

IF  Zl* 

=  "M"  THEN  368 

248 

IF  Zl* 

=  "G"  THEN  378 

260 

IF  Zl* 

a  "U"  THEN  948 

278 

IF  Zl* 

a  "L"  THEN  I860 

■*»  »  4L 

IF  Zl* 

a  "D"  THEN  CALL  -  151 

274 

IF  Zl* 

a  -C"  THEN  59 

U  (A  (A  t  j 
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280  HOME  :  GOTO  120 

90  INPUT  “WHAT  FIXED  VALUE  OF  U-PET?" *T1 
00  FO  =  409. 6s  T2  =  I NT  <T1  *  F0  ) : X  =  PEEK  <494 120 
20  FOR  TA  =  1  TO  T2:X  =  PEEK  < 49410):  NEXT  TA 
50  HOME  :  GOTO  12© 

360  X  =  PEEK  <49412):  HOME  :  GOTO  128 
378  INPUT  "FULL  UOLTAGE  RANGE? <V  OR  N)"SA* 

380  IF  A*  <  >  "V"  THEN  420 

390  U8  =  10: U9  =  0 :  GOTO  440 

420  INPUT  "MIN.  UOLTAGE  OF  RETARDER?" 5 U9 

438  INPUT  "MAX  UALUE?"SU0 

440  INPUT  "READINGS-- UOLTAGE  STEP?" 5 R 

450  INPUT  "HOW  MANV  LOOPS?" 5 L 

468  N  =  2.4414E  -  3*X  »  I NT  <<U0  x  N)  -  1):V  =  I NT  < 
U9  N ) :  G 1  »  X  -  VtVB  =  I  NT  <V  x  256) 

580  VB  =  I  NT  <V  x  256 ) :  VA  =  V  -  <256  *  VB):GB  =  I  NT 
CGI  x  256):  GA  =  G1  -  <256  *  SB) 

540  IF  81  >  4095  THEN  PRINT  “MAX  EXCEEDS  4095  COUNTS 

i  H 

545  IF  G1  >  4095  THEN  STOP 

550  POKE  27394, R:  POKE  27395, R:  POKE  27396, L:  POKE  27 
481, VA:  POKE  27402, VB:  POKE  27397, GA:  POKE  27398, 
GB 

620  PRINT  "THESE  UALUES  HAUE  BEEN  LOADED." 

630  CALL  27413 

640  MAX  =«  G1 :  HI  a  PEEK  <27409) 

660  S0  =  <  I NT  <<U1  -  MI)  *  409.6)  *  3  -*•  32739) 

670  IF  HI  =  1  THEN  PRINT  "":  PRINT  "♦ERROR  IN  DATA  B 

US  OR  AxD  POWER!!!!!*":  PRINT  “":  GOSUB  3000 

633  IF  HI  =  1  THEN  12© 

634  FOR  I  *  1  TO  5 

685  FOR  K  *  1  TO  8:  POKE  27634, K:  CALL  27685:  NEXT  K 

636  NEXT  I 

690  POKE  32763, R:  POKE  32769, L:  POKE  32770, VA:  POKE  3 
2771 , VB:  POKE  32772, GA:  POKE  32773, GB:  GOSUB  10 

691  INPUT  "U— TIP?" *  U1 

760  POKE  32774,  I  NT  <U1  -  I  NT  <U1  x  256)  *  256):  POKE 
32775,  I NT  <U1  X  256):  INPUT  "U- ANODE" 5 U2:  POKE  3 
2776,  I NT  <U2  -  I NT  <U2  x  256)  *  256): 

300  POKE  32777,  I NT  <U2  x  256) 

390  INPUT  "LENS  UOLTAGE" ; U5:  POKE  32783,  INT  <U5  x  2) 

:  INPUT  "TEMPERATURE  IN  KELUIN" S U7 
910  POKE  32734,  INT  <U7  x  4):  INPUT  "WHAT  OFFSET  <RET 
=0)?"?U3:  POKE  32736,  INT  <U3  x  256) 

923  POKE  32735,  INT  <U3  -  INT  <U3  x  256)  *  256): 

924  INPUT  “LIGHT  OR  DARK < 1  OR  0)"?A2 

925  GOTO  120 

940  INPUT  "WANT  FAULT  TOLERANCE?" 5 F* 

945  IF  F*  <  >  "V"  THEN  955 

950  INPUT  "WHAT  FAULT  TOLERANCE?" 5 IM 

955  INPUT  "WHAT  I?  "Si: 12  »  I  *  3:  FOR  II  *  12  TO  409 
6: PI  -  3  *  11 

960  IF  <32768  ♦  PI)  >  481 ?8  THEN  STOP 
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980  DEF  FN  UA<P2>  =  PEEK  <  P2  ♦  32789  +  Pt> 


PP 

=  FN  U A  <  0 

»  +  ( ^’56  * 

FN  ^A* 1 > > >  + 

<  256  '' 

< 

FN  UA<2>>> 

+  256  -  3 

* 

PEEK  <28679 

+  Il> 

991 

PQ 

=  FN  OA<3: 

•  256  * 

< 

FN  UA<4>  > >  + 

'•  256 

2  * 

r 

FN  UA<5>>> 

+  256  3 

* 

PEEK  <  28679 

+  n  + 

1  » 

993 

IF 

PQ  -  PP  < 

-  1 5000000 

THEN  POKE  C 

28679  + 

11  + 

1  ,  1 :  GOTO  990 

1001  PR  =  PQ  -  PPi  PRINT  <32788  +  P 1 > , PP , PR 

1002  IF  Ff  <  >  "V"  THEN  NEXT  II 

1003  IF  HBS  < PR >  >  IM  THEN  GOSUB  2040 

1005  PQ  =  FN  UA<3>  ■+•  <258  *  <  FN  UA  <  4  .'•>  >  >  +  <258  2  + 

<  FN  OH  < 5  > > > 

1010  PR  =  PQ  -  PPi  PRINT  <32788  +  p  1  } , PP , PR :  NEXT  11:  PRINT 
PRINT  "":  HOME  :  GOTO  128 

1080  PRINT  '*  ***.%***":  print  ’•  THIS  IS  A  DATA 

OIEWSrt  PLOT":  PRINT  "  OUER  ANY  RANGE  OF" 

1078  CLEAR 

1090  PRINT  "  THE  ENERGY  DISTRIBUTION.":  PRINT  " 
ft******":  PRINT  "":  GOSUB  1820 
1100  PL  =>  PEEK  <  32788  >  *  PEEK  <  32789  > 

1182  INPUT  "CLEAR  PREOIOUS  GRAPHS?" 5 Gl*:  IF  Gl*  ■  "V" 

THEN  HGR 

1185  INPUT  "HARDCOPY  ALSO?"?P*5  IF  P*  *  “V"  THEN  INPUT 
"READY?" 521*:  GOSUB  1908 

1188  POKE  -  18297.0:  POKE  -  18300-0:  POKE  -  18304, 

0:  HCOLOP=  3:  POKE  -  18301,0:  FOR  XI  *  0  TO  I NT 

<  N  > 

1178  X  =  I NT  <X1  *  N 1 > :  HPLOT  X, 157:  NEXT  XI 

1181  HPLOT  0,73  TO  273,78:  FOR  X2  =  0  TO  INT  <N  ✓  4> 

* X  =  INT  <X2  *  N1  *  4> :  HPLOT  X, 133:  HPLOT  X, 159 
:  NEXT  X2:  HPLOT  0,158 
1190  FOR  II  *  8  TO  278 

1202  2A  =  INT  <11  *  <ZC  -  T0>  ✓  27S> 

1203  T2  =  T0  +  2A  -  S3 

1204  T1  -  T2  +  2  *  S3 

1205  32  =■  S0  -►  3  *  ZA  -  3  *  S3 

1210  SI  -  S2  +  6  *  S3 

1211  PRINT  ”S1*"S1, "T1="T1 

1213  GOTO  1215 

1214  S2  -  S0  -  3  *  33:  T2  =■  T0  -  S3 

1215  DEF  FN  PP<S2>  »  PEEK  <S2>  258  *  PEEK  <S2  + 

1>  -*■  <258  A  2>  *  PEEK  <S2  2> 

1218  DEF  FN  PA<T2>  »  258  A  3  *  PEEK  <T2> 

1217  PP  *  0: PQ  *  0: A9  ■  PEEK  <6> 

1218  FOR  A  -  -  A9  TO  A9 

1220  PP  *  PP  +  FN  PP< 32  +  3  *  A>  +  FN  PA<T2  +  A> 

1230  PQ  ■  PQ  +  FN  PP<S1  3  *  A>  +  FN  PACT1  -t-  A> 

1 235  NEXT  A 

1240  PR  =  <PQ  -  PP>  ✓  <RL  *  AA> 

1241  IF  PR  >  M2  THEN  M2  »  PR:  POKE  32766, <  INT  <MZ  - 
256 > > *  POKE  32767, M2  -  256  *  <  INT  <MZ  --  256 >> 

1250  PS  =  PR  *  FI  /  S3 

1251  IF  23*  =  "V"  THEN  PS  »  PR  *  156  M2 
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1260 

IF 

PS  >  156 

THEN 

HPLOT  TO  I 1 

1261 

IF 

PS  >  156 

THEN 

PS  =  F'S  -  156 

1262 

IF 

PS  >  156 

THEN 

1261 

1270 

IF 

PS  <  0  THEN  P 

S  =  0 

1275 

IF 

P*  =  "V" 

THEN 

GOSUB  2000 

123®  PT  *  156  -  PS 

1290  HPLOT  TO  1 1 , PT 

1291  B5  =  B6 

1292  IF  F'S  <  73  THEN  B6  =  0 

1293  IF  PS  >  =  73  THEN  B6  =  1 

1294  IF  ABS  <B6  -  B5?  =  1  THEN  HPLOT  11,156  TO  11,1 
50:  HPLOT  1 1 , PT 

1300  NEXT  II 

1305  HI  =  0:  SO SUB  3000 

1310  INPUT  "GO  BACK  TO  MENU  CHIT  M?  OR  MORE  PLOTS'?"  ?  Z 
2$:  IF  Z2*  =  "M"  THEN  TEXT  :  HOME  :  CLEAR  :  GOTO 
120 

1330  GOTO  1060 

1350  PRINT  "R="  PEEK  < 32763 ?"  L="  PEEK  < 32769V  P*"<  PEEK 
<32763?  *  PEEK  <32769? > 

1355  US  =  PEEK  <32771?  *  256  +  PEEK  <32770? 

1360  UT  =  PEEK  <32774?  +  256  *  PEEK  <32775? 

1365  UA  =  PEEK  <32776?  +  256  *  PEEK  <32777? 

1370  PRINT  "U— START®" <  <  I NT  <US  /  4.096??  *  .01? 

1375  PRINT  "U— TIP* "UT"  U- ANODES "UA"  U-DIFF«”UT  - 

UA 

1330  PRINT  "OFFSET  UOLTAGE:"<  PEEK  <32735?  +  PEEK  <3 
2736?  *  256? 

1385  PRINT  "TEMPERATURE  IN  KELUIN: "4  *  PEEK  <32784? 

1390  PRINT  "LENS  UOLTAGE: "  PEEK  <32783?  *  2 
1450  PRINT  "DATE:"  PEEK  <32779?"/"  PEEK  <32780?"/"  PEEK 
<32731 ? 

1470  PRINT  "SAMPLE  TVPE:"  PEEK  <32773? 

1480  GET  KVs  GOTO  120 

1620  INPUT  "NORMALIZE?"? 23*: MZ  »  <256  *  PEEK  <32766? 

PEEK  <32767?? 

1626  IF  Z3*  =  "N"  THEN  MZ  *  0:  INPUT  "FACTOR  GUESS?"? 

FI 

1629  F0  =  409.6 

1630  MI  *  <  PEEK  <32771?  *  ?56  +  PEEK  <32770??  /  F© 

1632  MIX  »  1000  *  Ml: MI  a  MIX  /  1000 

1640  PRINT  "THE  DATA’S  STARTING  UOLTAGE  IS  "MI"  UOLTS 

M 

1645  INPUT  "HOW  MANV  POINTS  DO  VOU  WISH  TO  AUERAGE?" ? 

S3 

1647  INPUT  "ON  GRAPHING-WHAT  MESH  SIZE  FOR  AUERA6ING? 

" ?  A9*  POKE  6, A9 
1643  AA  *  2  *  a9  +  l 

1650  INPUT  "WHAT  UOLTAGE  DO  VOU  WISH  TO  START  AT?"?U1 

1659  T0  a  <  INT  <<U1  -  MI?  *  F0?  ♦  23679? 

1660  S0  a  <  INT  <<U1  -  MI?  *  F0?  *  3  >  32739? 

1661  INPUT  "WHAT  MAX  UOLTaGE" ?  ZB 


131 


1663  ZC  =  I NT  <  <2B  -  Ml?  *  F 8  >  +  28679 

1664  N  =  <ZB  -  Ul>  *  20!  HI  =  273  N 

1665  RETURN 

1961©  X  =  PEEK  <49438? :X  =  PEEK  <49439?!  FOR  I  =  1  TO 
100! X  =  PEEK  <49425?!  NEXT  I:  INPUT  "READY?" ? Z1 * 

1903  INPUT  "PLOT  T I CS? "‘Plfi  IF  Pit  *  "N"  THEN  RETURN 

1 904  X7  =  0 

1905  FOR  XI  =  0  TO  I NT  <N?:X  =  I NT  <H1  *  3  *  XI ?  - 
I NT  <  N 1  *  3  *  <X1  -  1?  ) 

1906  P9  =  10 

1907  IF  X7  =  0  THE '4  P9  =*  20 

1912  FOR  I  =  1  TO  200!  NEXT  I 

1914  FOR  I  =  1  TO  100:  NEXT  Is  FOR  I  =  1  TO  P9iV  =  PEEK 

<49428?!  NEXT  I 

1917  FOR  I  =  1  TO  100:  NEXT  Is  FOR  I  =  1  TO  P9s V  =  PEEK 

<49424?!  NEXT  I 

1925  FOR  I  =  1  TO  X!V  =  PEEK  <49426?!  NEXT  I 

1926  X7  =  X7  4!  IF  X7  -  16  *  0  THEN  X7  =  0 

1927  NEXT  XI 

1930  HI  =  0!  GGSUB  3000!  INPUT  "RAISE  PEN ! ! ! " 5 Z5* 

1931  X  *  PEEK  <49438? !X  a  REEK  <49439?!  FOR  I  =  1  TO 

1 00  s  X  =  PEEK  <49425?!  NEXT  I:  FOR  I  =  1  TO  10: X  = 
PEEK  <49424?:  NEXT  I 

1932  INPUT  " READY? ' >Zl*s  RETURN 

2000  A1  *  PS  -  Al:  IF  Al  >  0  THEN  A2  =  I  NT  <3  *  A1 

.5?:  IF  A2  <  >  0  THEN  FOR  I  =  1  TO  A2:X  =  PEEK 
<49424?:  NEXT  I 

2012  IF  Al  <  0  THEN  A2  =  INT  <  -  3  *  Al  +  .5?:  IF  A2 

<  >  0  THEN  FOR  I  =  1  TO  A2:X  =  PEEK  <49428?:  NEXT 

I 

2030  FOR  I  a  l  TO  3!X  =  PEEK  <49426?:  NEXT  I:A1  =  PS 
:  RETURN 

2040  POKE  <32768  +  PI  +  5?,<  PEEK  <32768  +  T1  2?? 

2041  POKE  <32768  -*■  PI  «•  4?,<  PEEK  <32768  +  PI  1?? 

2043  POKE  <32763  +  PI  +  3?,<  PEEK  <32768  +  PI?? 

2045  RETURN 

3000  FOR  I  »  l  TO  19:  IF  I  <  a  10  THEN  J  =  i 

3001  IF  I  >  10  THEN  J  =  20  -  I 

3002  A  «  is  POKE  27684, A:  CALL  27685!  NEXT  I 

3003  RETURN 


! 
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APPENDIX  B 

MACHINE  LANGUAGE  PROGRAM  FOR  DATA  ACQUISITION 

This  program  resides  in  memory  locations  *6B15  through  *6C88. 

It  is  written  for  the  Apple  II' a  Mostelc  6502  microprocessor.  The  left¬ 
most  column  of  the  printout  is  the  memory  location,  specified  in 
hexadecimal,  beginning  with  *6B15.  The  next  three  columns  hold  the 
actual  program  steps,  again  in  hexadecimal  notation.  The  rightmost  two 
columns  are  the  " disassembled"  mnemonic  notation  for  the  hexadecimal 
program  step  on  the  same  line. 

Lines  *6B16  through  *6B1C  cause  the  program  to  jump  to  several 
subroutines.  These  subroutines  set  initial  conditions  by  zeroing  all 
pertinent  data  storage  locations,  setting  the  retarder  to  its  minimum 
operating  voltage,  and  loading  the  R  and  L  counters  properly.  Steps 
*6B26  through  *6B3B  comprise  a  waiting  cycle,  and  test  the  ADC  for 
status.  If  conversion  is  not  complete  after  126  cycles,  control  is 
transferred  to  *6BA5,  which  sets  an  error  flag  and  exits  back  to  the 
Basic  program. 

Once  conversion  is  complete,  lines  *6B3D  through  *6B5F  perform 
a  quadruple  precision  addition  into  the  appropriate  memory  locations. 
Lines  *6B61  through  *6B42  increment  the  memory  locations  to  provide 
proper  sequencing  of  data,  and  test  the  R  and  L  loops  for  completion. 
Lines  *6C25  through  *6C52  generate  music  on  the  Apple' s  speaker. 

Lastly,  lines  *6C60  thorugh  *6C6F  test  the  data  bus  for  failure. 


r 
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*LL 


6B15- 

73 

SEI 

6816— 

20 

BO 

6B 

JSR 

<6BAD 

68 1 9— 

20 

04 

6C 

JSR 

<6C04 

6B1C- 

20 

09 

68 

JSR 

<6809 

6B1F- 

BO 

0C 

Cl 

LOB 

*C10C 

6B22— 

4C 

6*3 

SC 

JMP 

<6063 

6825- 

EB 

NOP 

6626- 

BD 

3B 

Cl 

LOB 

<C  1 0A 

6B29- 

B9 

79 

LOB 

#<79 

682B- 

30 

81 

6B 

STB 

<6801 

6B2E- 

CE 

01 

6B 

DEC 

<6801 

6B31- 

00 

03 

BNE 

<6836 

6B33- 

4C 

B5 

SB 

JMP 

<6BA3 

6B36- 

BO 

09 

Cl 

LOB 

<C- 1 09 

6B39- 

29 

F0 

AND 

#<F0 

6B3B- 

D0 

FI 

BNE 

<6B2E 

6B3D- 

B0 

00 

LDV 

#<00 

6B3F- 

BO 

0B 

Cl 

LOB 

<C10B 

6B42- 

BO 

03 

Cl 

LDB 

<C10S 

6B45- 

13 

CLC 

6B46- 

71 

F9 

ADC 

<<F9>, V 

6B43- 

91 

F9 

STB 

<<F9>,V 

6B4B- 

C3 

INV 

6B4B- 

BO 

09 

Cl 

LOB 

<C109 

6B4E- 

71 

F9 

ADC 

<<F9> . V 

6B50- 

91 

F9 

STB 

<<F9) > V 

6B52- 

C8 

INV 

6B53- 

B9 

00 

LOB 

#<00 

6B55- 

71 

F9 

BOC 

«F9>,V 

6B57- 

91 

F9 

STB 

<<F9>  »  V 

6B59- 

33 

DEV 

6B5B- 

38 

DEV 

6B5B- 

B9 

30 

LDB 

#<00 

6B3D- 

71 

ED 

BDC 

<<ED>*  V 

6B5F- 

91 

ED 

STB 

<<ed:>,  v 

6B61- 

CE 

83 

SB 

DEC 

<6683 

6B64- 

F0 

03 

BEQ 

<6869 

6B66- 

4C 

IF 

SB 

JMP 

<6B1F 

6B69- 

BD 

02 

SB 

LDB 

<6802 

6B6C- 

80 

03 

SB 

STB 

<6683 

Best  AvailzMe  toK 


I 
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6B6F— 

00 

02 

Cl 

LO0 

*C182 

6872— 

00 

12 

Cl 

LD0 

'■SCI  12 

6875- 

E6 

F9 

INC 

SF9 

6877- 

00 

02 

BNE 

S6B7B 

6879- 

E6 

F0 

INC 

SF0 

6878- 

E6 

F9 

INC 

SF9 

6870— 

00 

02 

BNE 

S688 1 

6B7F- 

E6 

F0 

INC 

SF0 

6881- 

E6 

F9 

INC 

SF9 

6883— 

00 

82 

BNE 

S6BS7 

6885- 

E6 

F0 

INC 

SFA 

6887- 

E6 

E0 

INC 

SEO 

6889- 

00 

02 

BNE 

S6830 

6888— 

E6 

EE 

INC 

SEE 

6880- 

C6 

EB 

DEC 

SEB 

6B8F- 

05 

E8 

LD0 

SE8 

6891- 

C9 

FF 

CMP 

#SFF 

6893- 

00 

01 

BNE 

S6866 

6895- 

C6 

EC 

DEC 

SEC 

6897- 

05 

EC 

LD0 

SEC 

6899- 

C9 

FF 

CMP 

#SFF 

6898- 

00 

C9 

BNE 

S6866 

6890- 

CE 

04 

6B 

DEC 

S6B04 

6800- 

F0 

06 

BEQ 

S6883 

6802— 

4C 

19 

68 

JMP 

S6B19 

6B05— 

EE 

11 

6B 

INC 

S6811 

6808- 

60 

RTS 

6B09- 

E0 

NOP 

6800- 

E0 

NOP 

6808- 

E0 

NOP 

6B0C- 

E0 

NOP 

6800- 

00 

00 

LDV 

#S00 

6B0F- 

82 

40 

LDX 

#S40 

6BB1- 

86 

FB 

STX 

SFB 

6BB3- 

82 

70 

LDX 

#S70 

6BB5- 

36 

F0 

STX 

SFA 

6BB7- 

82 

30 

LDX 

#S3© 

6B89- 

36 

FO 

STX 

SFO 

6BBB- 

02 

15 

LDX 

#S15 

6BBD- 

36 

FC 

STX 

SFC 
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6BBF- 

34 

F9 

STV 

♦F9 

6E'C 1  — 

98 

TVA 

6BC2- 

30 

1  1 

6B 

STA 

♦6B1  1 

6BC5— 

91 

F9 

3TA 

<SF9':> 

6BC7— 

E6 

F9 

INC 

♦F9 

6BC9- 

00 

FA 

BNE 

♦6BC5 

6BCB— 

E6 

FA 

INC 

♦FA 

6BCD— 

C6 

FB 

DEC 

♦FB 

6BCF- 

00 

F4 

BNE 

♦6BC5 

6BD1- 

A9 

40 

LOA 

4*40 

6BD3- 

35 

FB 

STA 

♦FB 

6805— 

60 

RTS 

6BD6— 

EA 

NOP 

6607- 

EA 

NOP 

6BD3- 

EA 

NOP 

6BD9- 

AO 

0A 

6B 

LOA 

♦6B0A 

6BDC- 

35 

FF 

STA 

♦FF 

66DE- 

AO 

09 

6B 

LOA 

♦6809 

6BE1- 

35 

FE 

STA 

♦FE 

6BE3- 

AD 

04 

Cl 

LOA 

♦Cl  04 

66E6- 

AD 

IE 

Cl 

LDA 

♦CUE 

6BE9- 

AO 

02 

Cl 

LOA 

♦Cl  02 

6BEC- 

AD 

12 

Cl 

LDA 

♦Cl  12 

6BEF- 

C6 

FE 

DEC 

♦FE 

6BF1- 

A5 

FE 

LDA 

♦FE 

6BF3- 

C9 

FF 

CMP 

#*FF 

6BF5- 

00 

F2 

BNE 

♦6BE9 

6BF7- 

C6 

FF 

DEC 

♦FF 

6BF9- 

A5 

FF 

LDA 

♦FF 

6BFB- 

C9 

FF 

CMP 

#*FF 

6BFD- 

00 

EA 

BNE 

♦6BE9 

6BFF- 

60 

RTS 

6C00- 

EA 

NOP 

6C01- 

EA 

NOP 

6C02- 

EA 

NOP 

6C03- 

EA 

NOP 

6C04- 

A5 

FC 

LDA 

♦FC 

6C06- 

35 

F9 

STA 

♦F9 

6C08- 

A5 

FD 

LDA 

♦FD 

6C0A- 

35 

FA 

STA 

♦FA 

Best  Available  Copy 


6000- 

39 

40 

LD3 

4*40 

6C0E- 

35 

13 

ST  3 

*13 

60 1 0- 

30 

05 

6B 

LD3 

*6605 

6013- 

35 

EB 

ST3 

*EB 

6015- 

3D 

86 

6B 

LD3 

*6B06 

60 18¬ 

35 

EC 

ST3 

SEC 

60 18¬ 

39 

70 

LD3 

4*70 

60 10- 

35 

EE 

ST3 

*EE 

60 1E- 

39 

07 

LD3 

4*07 

6020 — 

35 

ED 

ST3 

*ED 

bC2^,'- 

60 

RTS 

6023- 

E3 

NOP 

6024- 

08 

PHP 

6025- 

32 

FF 

LDX 

4*FF 

6027- 

3D 

30 

00 

LD3 

*0830 

6023- 

3D 

24 

60 

LD3 

*6024 

6C2D- 

20 

33 

FC 

JSR 

*FC3S 

6030- 

03 

DEX 

6031- 

D0 

01 

BNE 

*6034 

6033- 

60 

RTS 

6034- 

40 

27 

60 

JMP 

*6027 

6037- 

60 

RTS 

6033- 

10 

10 

BPL 

*6043 

6033- 

09 

18 

0R3 

4*18 

6030- 

09 

03 

0R3 

4*83 

6C3E- 

09 
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